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Abstract. We describe a method of disentangling the composite, 0.12-5 j-im continuum of symbiotic binaries. The observed 
SED is determined by the IUE/HST archival spectra and flux-points corresponding to the optical UBVRI and infrared JHKLM 
photometric measurements. The modeled SED is given by superposition of fluxes from the cool giant, hot stellar source and 
■ nebula including the effect of the Rayleigh scattering process and considering influence of the iron curtain absorptions. We 

. applied this method to 21 S-type symbiotic stars during quiescence, activity and eclipses. We isolated four main components 

£N| ■ of radiation and determined their properties, (i) Stellar radiation from the giant corresponds to a unique luminosity class - 

normal giants. Characteristic luminosities are 1 600 + 200 and 290 ± 30 Lq for red and yellow giants, respectively in our sample 
of objects, (ii) Hot object radiation during quiescence consists of the nebular and stellar component. The former radiates at a 
mean electron temperature of 19 000 K and its amount of emission suggests a mass-loss rate from giants via the wind at Mw = a 
few x 10~ 7 Mq yr~' . Radiation of the latter conforms well with that of a black-body photosphere at a characteristic temperature 
of 105 000 K. The corresponding effective radii are a factor of ~ 10 larger than those of white dwarfs, which thus precludes 
observing the accretor's surface. Extreme cases of AX Per and V443 Her, for which the hot star temperature from the fit is not 
capable of producing the nebular emission, signal a disk-like structure of the hot stellar source even during quiescence, (iii) Hot 
object radiation during activity consists of three components - the stellar and the low- and high-temperature nebular radiation. 
The stellar radiation satisfies that of a black-body photosphere at a low characteristic temperature of ~ 22 000 K (we call it 
the 1st type of outbursts) or at a very high characteristic temperature of a 165 000 K (2nd type of outbursts). All the active 
objects with a high orbital inclination show features of the lst-type of outbursts (here Z And, AE Ara, CD-43° 14304, TXCVn, 
BFCyg, CHCyg, CICyg, ARPav, AX Per), while AGDra represents the 2nd-type. The presence of a two-temperature type 
of UV spectrum and an enlargement of effective radii of the stellar source by a factor of ~ 10 with respect to the quiescent 
values during the lst-type of outburst suggest an expansion of an optically thick medium at the orbital plane in the form of 
a disk. The low-temperature nebula radiates at a mean electron temperature of 14 000 K and is subject to eclipses, while the 
high-temperature nebula, which is seen during eclipses as the only component, is characterized by T e > 30 000 K. 

Radiative and geometric properties of the main sources of radiation allowed us to reconstruct a basic structure of the hot 
object during the lst-type of outburst. There is an edge-on disk around the accretor. Its outer flared rim represents a warm 
pseudophotosphere of the hot stellar source, whose radiation is Rayleigh attenuated and affected by the iron curtain absorptions 
in the neutral gas concentrated at the orbital plane. The low-temperature nebula is placed just above/below the disk with a 
concentration at its edge as to be subject to eclipses and to 'see' well the central ionizing source. High above/below the orbital 
plane, there is a hot nebular emitting region. 

Key words, method: data analysis - stars: binaries - stars: symbiotics - accretion: accretion disks 



1. Introduction 

Symbiotic stars are interacting binary systems consisting of a 
cool giant and a hot compact star. Typical orbital periods are 
between 1 and 3 years, but can be significantly larger. The mass 
loss from the giant represents the primary condition for inter- 
action between the binary components. A part of the material 
lost by the giant is transferred to the more compact compan- 
ion via accretion from the stellar wind or Roche-lobe overflow. 
This process generates a very hot (7h ~ 10 5 K) and luminous 
(Lh ~ 10 2 - 1Q 4 Lq) source of radiation. How the generated 
energy is liberated determines two stages of a symbiotic bi- 



nary. During quiescent phases the hot component releases its 
energy approximately at a constant rate and spectral distribu- 
tion. The hot radiation ionizes a fraction of the neutral giant's 
wind, which gives rise to nebular emission comprising numer- 
ous lines of high excitation/ionization and the continuum. As 
a result the spectrum consists of three components of radia- 
tion - two stellar and one nebular. For eclipsing systems, the 
hot star radiation can be attenuated by Rayleigh scattering at 
the position where the mass losing giant is in front (Isliker et 
al. 1989). During active phases the hot component radiation 
changes significantly, at least in its spectral distribution, which 
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leads to the 2-3 mag brightening of the object in the optical. A 
common feature of active phases is a high-velocity mass out- 
flow, which can significantly disrupt the ionization structure of 
symbiotic nebulae. In some cases the nebular emission disap- 
pears entirely at the optical maximum (e.g. BFCyg), in others 
dominates the UV/optical spectrum (e.g. AGDra). 

Many particular aspects of this general view have been re- 
cently discussed and summarized with outlined problems by 
Corradi et al. (2003). Current pivotal problems in the sym- 
biotic stars research are connected mainly with understanding 
the nature of outbursts and the observed accompanied phenom- 
ena (Sokoloski et al. 2002; Bode 2003; Skopal 2003c; Lee 
2000). In this respect the nature of hot objects in symbiotic 
stars play a crucial role. A way to determine their fundamen- 
tal parameters is by modeling the ultraviolet spectrum. A pi- 
oneering work in this direction has been made by Kenyon & 
Webbink (1984,hereafter- KW), who elaborated a grid of syn- 
thetic spectra for various binary models. By comparing the cal- 
culated and observed ultraviolet continuum colours they sug- 
gested the type of the binary model (MS/WD accretor or a sin- 
gle hot star) with the corresponding accretion rate and the hot 
star temperature. Miirset at al. (1991, hereafter - MNSV) re- 
fined their method by developing a more sophisticated model 
of the nebular continuum accounting for fractional rates of pho- 
tons capable of ionizing He + and H. This approach requires in- 
cluding a geometrical structure of He ++ and H + ionized zones 
in the binary. They fitted the observed IUE spectra scaling the 
hot star contribution, its temperature (calculated as the Zanstra 
temperature) and a parameter X, which defines the geometry of 
the nebula in the sense of the STB ionization model (Seaquist 
et al. 1984, hereafter - STB). Both groups assumed the elec- 
tron temperature to be constant throughout the nebula, but kept 
its quantity at 10 000 K and/or 20 000 K (KW) or approximated 
it by values in a rather narrow interval of 12 800 - 15 000 K for 
the H + region (MNSV). 

The aim of this paper is to introduce a simple method of dis- 
entangling the composite SED in the UV/optical/IR continuum 
of symbiotic binaries (Sect. 3). In our approach we fit observa- 
tions by a model with 6 basic free parameters - temperatures of 
the three main emitting sources and scalings of their contribu- 
tions. In addition we consider effects modifying the continuum 
- the Rayleigh scattering process and veiling the continuum 
by a forest of blended Fe n absorptions. In the second part of 
this paper (Sect. 4) we apply the method to S-type symbiotic 
binaries during quiescence, activity and eclipses. Particularly 
we discuss cases in which the modeled parameters lead to a 
conflict with the observed properties. In Sect. 5 we summarize 
common characteristics of the isolated components of radiation 
and their sources. 



tra taken by the FOS on the board of the HST satellite. We 
selected spectra with well exposed continuum, taken with the 
large aperture and the short- and long-wavelength prime dur- 
ing one shift. We used the spectra taken during quiescence 
(if available, at dates of different conjunctions of the binary 
components), activity and optical eclipses. Main sources of the 
data covering the infrared region are represented by broad-band 
multicolour photometric observations published by Munari et 
al. (1992),Kamath& Ashok (1999), Taranova (2000), Glass 
& Webster (1973) and Feast et al. (1977). Measurements in 
the Ic and Rc bands of the Cousins system were converted into 
the Johnson system according to Bessell (1983). In addition, 
we summarized appropriate parts of light curves (LC) to dis- 
tinguish stages of quiescence and activity at the time of spec- 
troscopic observations. The main sources of the data here are 
represented by UBVR photometric observations collected dur- 
ing the last 15 years of our campaign for long-term monitoring 
of symbiotic stars (Skopal et al. 2004, and references therein). 
Other data for individual objects are referred in relevant sec- 
tions. Stellar magnitudes were converted to fluxes according to 
the calibration of Lena et al. (1999). We dereddened observa- 
tions with appropriate value of £b-v using the extinction curve 
ofCardellietal. (1989). 

In Table 1 we list the objects, ephemerides, colour excesses 
and distances we use in this paper. 

3. Disentangling the composite spectrum 

3.1. Model of the continuum 

According to the three-component model of radiation (Sect. 1) 
we can express the observed continuum flux, F(A), of symbi- 
otic stars as a superposition of two stellar components, Fh(A) 
and F g (A), from the hot and cool giant star, respectively, and 
the nebular contribution, F N (A), as 

F(A) = F g (A) + F h (A) + F N (A) (1) 

in units of erg crrT 2 s _1 A -1 . We approximate radiation from 
the giant by an appropriate synthetic spectrum, F*! ynth '(7eff)> 
scaled to flux points given by the broad-band photometry. Thus 
we can write 



FJA) = Ff ath XTcs), 



(2) 



The effective temperature of the giant, T t s, and the scaling fac- 
tor of its spectrum determine the bolometric flux at the Earth's 
surface as 

F° bs = jF^ th (T eB ) dA = k g JnB A (T eB ) dA = k g crT 4 eff (3) 



2. Observations 

We apply our method (Sect. 3) to well observed S-type symbi- 
otic stars, for which the low-resolution ultraviolet spectroscopy 
and the broad-band optical and infrared UBVRI and JHKLM 
photometry are available. The most important data source we 
used to cover the ultraviolet domain is the Final archive of the 
IUE satellite. In the case of AG Peg we analyzed also the spec- 



in units of erg cm s .In terms of the giant's luminosity L g = 

S u 1 eff ' 

U iR e \ 2 



4nRgCrT^ s , the scaling factor, k g , can be expressed as 



(?)■* 



(4) 



where 8 g is the angular radius of the giant. We approximate the 
hot star continuum by a black-body radiation at a temperature 
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r h . In addition, we consider that it can be attenuated by neu- 
tral atoms of hydrogen due to the Rayleigh scattering process. 
Thus, we express the second term on the right side of Eq. (1) 
as 



F h (A)=k h xnB /l (T h )e- n ^, 



where = F^AJ/nB^T^) is a dimensionless scaling factor, 
which can be expressed in the same form as that for the giant's 
radiation, i.e. 



4nd 2 crT* 



t) ^ 



(6) 



where Lh, 7h, and $h are the luminosity, temperature, ef- 
fective and angular radius of the hot stellar source, respec- 
tively. Further in Eq. (5), «h [cm -2 ] is the column density of 
H atoms and <x^ [cm 2 ] is the Rayleigh scattering cross-section 
for atomic hydrogen. According to Nussbaumer et al. (1989) 
we calculate the scattering cross-sections as 



z 



fu 



(A/A lk ) 2 - 1 



(7) 



where cr e = 6.65E-25 crrT 2 is the Thomson cross-section, f\ k 
are oscillator strengths of the hydrogen Lyman lines and Au are 
the corresponding wavelengths. An interesting feature (from 
the astrophysical point of view) of the cr^ function is its rapid 
decrease from Lya to ,11055 A (see Fig. 2 of Nussbaumer et 
al. 1989). This effect should be seen as a rapid increase of 
the far-UV continuum beyond the Lya line to its scattering free 
values at ,11055 A. For example, the SWP spectra from the IUE 
archive begin at ~ 1 140 A, which means that the effect could be 
detected only for rather low values of «h ~ 10 23 cm 2 and on 
the well exposed spectra. 

The last term on the right side of Eq. (1) represents flux in 
the continuum emitted by the symbiotic nebula. For the purpose 
of this paper we consider nebular contributions from hydrogen 
and doubly ionized helium generated by processes of recom- 
bination (f-b transitions) and thermal bremsstrahlung (f-f tran- 
sitions). We neglect contributions from singly ionized helium, 
because its emission coefficient is comparable, in both the pro- 
file and the quantity, with that of hydrogen, but it is much less 
abundant in the ionized medium (to about 10%). This similar- 
ity also makes it difficult to identify any signatures of the emis- 
sion from the singly ionized helium in the observed spectrum. 
Therefore, and also for simplicity, we express the nebular flux 
as 



F N (A) 



\nd 2 \ Y 



e A (H, T e )N n + + 



s A (Ue + , T e )AW + N e dV, 



(8) 



where d is the distance to the object and s A (erg cm s A ) is 
the volume emission coefficient per electron and per ion under 
consideration, which includes all acts of f-b and f-f transitions. 
Nu+, A^He ++ and N s are concentrations of the ionized hydro- 
gen (protons), doubly ionized helium and the electrons, respec- 
tively, and r e is the electron temperature. Further we use the 



following two assumptions for the symbiotic nebula: (i) Both 
the H + and the He ++ regions occupy the same volume (in con- 
trast to the usual case, for which the volume Vh+ > VW+)- 
Then the ions in such nebula are in the pool of electrons with 



(5) N e = (l+2a)Nw 



where 



a = AW+/N H + 



(9) 



(10) 



represents an average abundance of the doubly ionized helium 
throughout the nebula, (ii) The nebular radiation is character- 
ized by a uniform electron temperature, L e , which is constant 
within the nebula. The assumption (i) is based on that the spec- 
tra include contributions from all the present ions irrespective 
of their displacement within the nebula. Then the average abun- 
dance a represents a lower limit of that in the real He ++ re- 
gion. Specially, in the case that Vh+ = VW+, also a = a(He ++ ). 
Contrarily, if Vn+ » VW+ then a <K a(He ++ ), which implies 
that the contribution from the He ++ region is very small with 
respect to that from H + region. These limiting cases can be 
distinguished from observations. For example, a large quantity 
of a derived from observations (to say around 0.1, which can 
be considered as a maximum given by the helium abundance 
fl(He)) means that the helium is doubly ionized everywhere in 
the nebula, i.e. VHe++ ~ Vh* ■ The assumption (ii) is verified by 
good fits of the continuum at spectral regions, where the He ++ 
and H + nebular continuum dominate the spectrum (see Figs. 2 
- 22). So, with the aid of Eqs. (9, 10) and the assumption of 
uniform L e in the nebula, the expression (8) takes the form 



EM 

F n(A) = s a (R, He + , T e , a), 



(11) 



where the emission measure EM = J y A^ + dV and the emission 
coefficient 



«i(H, He + , T e ,a) = (1 + 2a)e A (U, T e ) + 

2(1 + 25)^(He + , T e ). 



(12) 



For small values of a, to say a < 0.05, s A =e A (¥L, T e ) + 
fl£i(He + , r e ), which can be applied, for example, to symbi- 
otic nebulae during quiescent phases that are characterized by 
> ^He++- Finally, according to relations (2), (5) and (11), 
we can express Eq. (1) in the form 



Fx = F^\T tS ) + 



k h XnB A (T h )e- n ^ + 
k N xs A (H,He + ,T e ,a), 



(13) 



where £n is the scaling factor for the nebular component of 
radiation. According to Eq. (11) 



k N 



EM 
4^ 



[cm" 5 ] 



(14) 



3.1 .1 . Influence of the line blanketing effect 

This effect represents an additional source of the absorbing 
spectrum due mostly to Fen transitions. This so-called "iron 
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SWP38698+LWP 17832 (29/04/90) 
MWL 560 SWP38361+LWP17534 (14/03/90) 

SWP42580+LWP21366 (29/09/91) 
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Fig. 1. IUE spectra of MWC 560 obtained on 29/04/90, 14/03/90 and 
29/09/91 with different amounts of absorption by the veiling gas. 



curtain" was first identified by Shore & Aufdenberg (1993) as 
the origin of the anomalies in the emission-line fluxes result- 
ing from the differential absorption by the environment as the 
path length along the line-of-sight changes. Generally, the in- 
tervening absorption originates in a slab of cool gas that veils 
the source of the hot radiation. In symbiotic binaries the veil- 
ing gas can be represented by the neutral wind of the giant, the 
effect of which is largest at/around the giant's inferior conjunc- 
tion (e.g. Dumm et al. 1999). In this study we find that such a 
slab of absorbing gas can temporarily be created at the orbital 
plane during the Z And type of outbursts. 

The line blanketing effect results in a complex profile of 
the observed continuum. On the low-resolution spectra broad 
absorption bands at 1500 to 1800 A and 2 300 to 2 800 A ac- 
companied by spurious emissions are characteristic features. 
To estimate the level of the continuum at wavelengths with the 
lowest absorption, we have chosen three representative low- 
resolution IUE spectra of MWC 560 to obtain a template for a 
low, modest and strong influence of the continuum by the iron 
curtain (Fig. 1). With the help of these observations and con- 
sidering theoretical calculations given by Shore & Aufdenberg 
(1993), Home et al. (1994) and Dumm et al. (1999) we se- 
lect points of the continuum at ~ 1280, ~ 1340, ~ 1450, ~ 1500, 
(-1580-1600), around 1800, - 1950-2150, (-2450), -2640 
and beyond 2800 A depending on the amount of absorption. 

The resulting effect of the absorbing gas in the spectrum 
of symbiotic binaries depends also on its relative position with 
respect to the location of different emitting regions in the sys- 
tem. Generally, radiation from the hot stellar source is absorbed 
most, while the nebular emission - having its origin in a region 
located well away from the hot star - is less affected. An addi- 
tional complication is the superposition of the strong emission 
lines on the broad absorption bands of the iron curtain. For ex- 
ample, emission of Hen 1640, Oi 1641] and Ora 1664 A lines 
fills up in part the underlying absorption bands, which results in 
an artificial absorption located at/around 1700 A. Its depth re- 
flects the magnitude of the iron curtain effect. This part of the 
spectrum cannot be used to estimate the level of the continuum. 



3.2. The model 

To determine the giant's continuum we scaled a synthetic spec- 
trum to the broad-band [BVRI]JKLM photometric fluxes (de- 
pending on the giant's contribution in the optical). This ap- 
proach is applicable for S-type symbiotic stars (i.e. with a stel- 
lar type of the spectrum; no significant dust emission is present 
in the IR wavelengths). An appropriate synthetic spectrum, 
which best matches the photometric flux-points, was selected 
from a grid of models made by Hauschildt et al. (1999). We 
found it satisfactory to use only a (r eff , log(g)) grid and to 
smooth the spectrum within 20 A bins to keep a basic profile 
of the giant's continuum. The resulting parameters are T e tf and 
g , which define the observed bolometric flux, F° bs (Sect. 3.1). 
For most of the objects we selected log(g) = 0.5, for CH Cyg 
and T CrB we put log(g) = and for yellow symbiotics, 1.5- 
2, in accordance with their detailed atmospheric analysis (e.g. 
Smith etal. 1997). 

Having determined the radiation from the giant we sub- 
tracted its contribution from the selected flux-points of the UV- 
continuum, F \, s (Ai), and fitted them by the synthetic continuum 
expressed by Eq. (13), for which the function 

1 N 

i=l 

reaches a minimum. In this way we determined the other 6 free 
parameters, k^, Th, riu, k N , T e and a, which define the individual 
components of radiation. In Eq. (15), F(Ai) are the theoretical 
fluxes given by Eq. (13) at selected wavelengths A[. N is their 
number. We estimated fluxes F b s (Ai) at about 10 to 20 wave- 
lengths between 1 250 and 3 300 A depending on quality of 
the used spectrum. Because of noise, numerous emission lines 
and absorption features we estimated the continuum fluxes by 
eye at the wavelengths specified in section 3.1.1. Then, if pos- 
sible and applicable, we subtracted the contribution from the 
hot temperature nebula (see Sect. 5.3.1). Second, we prepared 
a grid of models for reasonable ranges of the fitting parame- 
ters and selected that corresponding to a minimum of the x 2 
function. We found the initial ranges for k h , T h and n H param- 
eters by comparing the Planck function for different temper- 
atures (15 000 - 40 000 K for some active phases) attenuated 
with n H ~ 10 23 ctrT 2 (if applicable) to the far-UV spectrum 
between about 1 250 and 1 600 A. For the spectra from quies- 
cent phases we adopted in some cases their Zanstra tempera- 
ture according to MNSV and/or determined their lower limits 
as introduced in Sect. 4.1. Section 5.2.2 gives more details. The 
range for T e can be estimated according to the slope of the near- 
UV continuum from about 2400 A to the Balmer jump, where 
the nebular emission represents a dominant contribution. The 
average abundance of the doubly ionized helium, a, can be es- 
timated from the difference of the continuum level in the range 
of about 1600 A to 1980 A and about 2 300 A to -3 200 A. 
If we use the IUE spectra, these parts are also well defined. 
Unfortunately, the part from the beginning of the LWP(R) spec- 
tra to about 2 300 A is very noisy and the following region in 
between - 2 300 and - 2 700 A is often affected by the iron 
curtain absorptions, which makes it difficult to recognize sig- 
natures from the He ++ emission. Finally, the scaling factor fc N 
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results from the fitting procedure. If the initial grid of models 
(one model = one combination of the fitting parameters) in- 
cluded the minimum of the^ 2 function, then a close grid could 
be determined around the minimum. Repeating this approach 
we selected models fitting the observed fluxes within their un- 
certainties. 

3.3. Uncertainties of the fitting parameters 

To estimate uncertainties of the resulting fitting parameters, we 
determined axmux quantity, which separates the reliable models 
from those given by all combinations of the grid parameters. 
The models with^f < Xmax then reconstruct the observed fluxes 
within their uncertainties. According to evaluation of the well 
exposed spectra from the IUE archive, they are of about 10% 
of the measured values. By preparing many trials we found an 
approximative relation for^ max as 

AWx ~ 0.05 x F obs , (16) 

where F D \, S is the average of the F b S (/li) fluxes. The Xmax 
value then bounds the ranges for the fitting parameters in the 
X(nH, Th, kh, 5, T e , fc N ) diagrams. Note that a change of the pa- 
rameter under consideration within the range of its possible 
values requires a different combination of the other fitting pa- 
rameters to keep the x value of the corresponding model un- 
der the Xmux limit. The primary source of uncertainties is the 
accuracy of the measured continuum. The ranges of possible 
values of the fitting parameters are also given by dependen- 
cies of the synthetic continuum on individual parameters. The 
profile of the far-UV continuum depends strongly on the n H pa- 
rameter (Eq. 5). Any other parameter can exceed its influence 
on the continuum at these wavelengths. For the cases with a 
significantly Rayleigh-attenuated far-UV continuum, the^-(n H ) 
function has a rather sharp minimum. A complication here is a 
lower temperature of the stellar source, 7j, ~ 20 000 - 30 000 K 
measured during some active phases, because the maximum 
of the corresponding Planck's function lies around the far-UV 
region or just beyond it at shorter wavelengths. In this case a 
higher temperature requires also a slightly higher « H to keep 
the fit under Xmux- This dependence mainly enlarges the range 
of possible values of T b . The profile of the nebular continuum 
from about 2 000 A to the long-wavelength end of the spec- 
trum depends significantly on the parameters a and T e . For the 
spectra strongly affected by the iron curtain absorption it is not 
possible to determine unambiguously the a parameter. In addi- 
tion, our models (Sects. 4.2 - 4.22) suggest it to be very small. 
In conclusion, ranges of the fitting parameters Th, «h> T e , kh 
and k N were estimated to be 20% to 50% of their best values. 

Uncertainties for the stellar component of radiation from 
the cool giant are given by those in T e g and the scaling. The 
former is proximately +100K (see also Skopal 2000) and the 
latter is given by that of photometric measurements, which we 
assumed to be less than 10%. Then the relative mean square 
error in the observed bolometric flux from the giant is about of 
15%. 



Table 1. Ephemerides, reddening and distances 



Object 


Sp. conj." 


^orb 


£b-v 


d 


Ref. 




[MJD] 


[day] 


[mag] 


[kpc] 




EG And 


50683.2 


482.57 


0.05 


0.59" 


1,2 


ZAnd 


14625.2 


757.5 


0.30 


1.5 


3,4 


AE Ara 


50217 


812 


0.25* 


3.5* 


5 


CD-43 C 


45 567 


1448 


0.2 


2.1* 


6 


TCrB 


47 861.7 


227.57 


0.15 


0.96 


1,7,8 


TXCVn 


45 130.5 


198 


0.0 


1.0* 


9 


BFCyg 


45 395.1 


757.2 


0.35 


3.8* 


10 


CHCyg 


45 888 


756 


0.0 


0.27 


11,12 


CICyg 


11902 


855.25 


0.35 


2.0* 


13,14 


V1329Cyg 


27 687 


958.0 rf 


0.37 


4.2* 


15 


LTDel 


45 910 


478.5 


0.20* 


3.9* 


31 


AGDra 


43 629.2 


549.73 


0.08 


1.1* 


16,17 


CQDra 


42932 


1703 


0.10 


0.178 


18,19 


V443 Her 


43 660 e 


594 


0.25* 


2.2* 


20 


YY Her 


40637 


592.8 


0.20 


6.3* 


3,21 


RWHya 


49 512 


370.4 


0.10 


0.82* 


22,2 


SYMus 


50176 


625.0 


0.35* 


1.0* 


23 


ARPav 


11266.1' 


604.45 


0.26 


4.9 


24-6 


AG Peg 


27 664.2 


812.6 


0.10 


0.80 


3,27 


AX Per 


36 673.3 


679.9 


0.27 


1.73 


28-30 


IVVirS 


49 016.9 


281.6 


0.20* 


1.3* 


32 



a - inferior conjunction of the cool component 

b - this paper 

c-CD-43° 14304 

d - pre-outburst ephemeris 

e — J £>Min(U) 

/ - a mean linear ephemeris of eclipses 
£-BD-21°3873 

References: 1 - Fekel et al. (2000a), 2 - MNSV, 3 - Skopal (1998), 
4 - Mikolajewska & Kenyon (1996), 5 - Mikolajewska et al. (2003), 
6 - Schmid & Nussbaumer (1993), 7 - Selvelli et al. (1992), 8 - 
Belczynski & Mikolajewska (1998), 9 - Kenyon & Garcia (1989), 
10 - Fekel et al. (2001), 11 - Hinkle et al. (1993), 12 - Viotti et al. 
(1997), 13 -Aller (1954), 14 - Kenyon et al. (1991), 15 - Schild & 
Schmid (1997), 16 - Galis et al. (1999), 17 - Birriel et al. (2000), 18 
- Reimers et al. (1988), 19 - Perryman et al. (1997), 20 - Kolotilov 
etal. (1995), 21 - Munari et al. (1997a), 22 - Schild et al. (1996), 
23 - Dumm et al. (1999), 24 - Skopal et al. (2001b), 25 - KW, 26 - 
Schild etal. (2001), 27 - Kenyon et al. (1993), 28 - Skopal (1991), 
29 - Mikolajewska & Kenyon (1992a), 30 - Skopal (2000), 31 - 
Archipova et al. (1995), 32 - Smith et al. (1997) 

4. Application to S-type symbiotic stars 

4.1. Fitting and derived parameters 

Fitting parameters, k g , kh, k^ and corresponding temperatures 
determine stellar luminosities and emission measure of the neb- 
ula (Eqs. 4, 6, 14). A useful parameter is the ratio of the hot star 
and the cool giant luminosity, which is independent of the dis- 
tance. According to our notation this can be expressed as 

^ = (17) 
L g k g \ r e ff / 

At the wavelengths, where the nebular emission dominates the 
spectrum, one can estimate the upper limit of EM directly from 
observations by assuming that F n(/1) = F bs(X), which is usu- 
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ally fulfilled for fluxes in the near-UV region. Then, according 
to Eq. (1 1) the emission measure is limited by 



EM = And 



2 F obs (A) 



s A (U, He + ,r e ,a) 
We calculate the luminosity of the nebula from Eq. (1 1) as 



(18) 



S+OO 

Ln = EM ( ei(H, He + , T e , a)dA, 
J912 



(19) 



which is valid for an ionised medium that is optically thick in 
the Lyman continuum. Finally, we introduce a parameter 



<5 = 



EM obs 
EM B ' 



(20) 



where EM obs is the emission measure derived from the model 
(Eq. 14) and EM B = Lph/<*B results from the equilibrium equa- 
tion between the rate of ionizing photons, L pb (photons s -1 ), 
and the rate of ionization/recombination acts within the ion- 
ized region; as (cm 3 s _1 ) stands for the total hydrogen recom- 
bination coefficient (STB). Thus EMg represents a maximum 
of the nebular emission, which can be produced by the ioniz- 
ing source under conditions of case B in the hydrogen plasma. 
This implies that 8 < 1 . Values close to 1 indicate that all ion- 
izing photons are converted by the circumbinary medium and 
are not lost. Small values of 8 indicate the opposite case. Using 
the relation for L pb from Skopal (2001a) the 8 parameter can 
be expressed as 



8=% B /f(T b ), 



where the function 



- r\ 

he Jo 



AB A (T h )dA. 



(21) 



(22) 



For example, /(8 10 4 K) = 4.6E+25, f(l 10 5 K) = 1.1E+26, 
/(1.3 10 5 K) = 2.7E+26and/(1.5 10 5 K) = 4.3E+26 cnrV 1 . 
The limiting case, 8 = 1, determines a lower limit of the 
temperature of the ionizing source, T™ n , required to produce 
the observed emission measure. Thus, the number of ioniz- 
ing photons produced by the hot object radiating at r™ n just 
balances the number of recombinations in the continuum. In 
modeling the ultraviolet continuum by Eq. (13) we often used 
Th = Tzanstra (from MNSV) for a first-model-fit. If the ionizing 
capacity of the hot stellar source was not consistent with the 
observed nebular emission (8 > 1) we searched for T™ n . As a 
temperature change in the range of T b > 10 5 K practically does 
not influences the resulting fit, we can get T h = r™ n by solving 
equation 



-a B -f(T h ) = 



(23) 



k b (T h ) 

for parameters fe N and T e from the first-model-fit. 

In the following sections we give results for individual 
objects with special attention to common properties of their 
SEDs during quiescent and active phases as well as eclipses. 
The resulting parameters of the reconstructed continuum are in 
Tables 2, 3 and 4. Corresponding synthetic continua together 
with observations are shown in Figs. 2 - 22. 



4.2. EG Andromedae 

EG And is a quiet symbiotic star - no active phase has been 
recorded to date. UBV LCs display a periodic orbitally-related 
variation of a double-wave in the profile (Fig. 2). Occasionally, 
0.3 - 0.5 mag flares are observed in the U-LC. The infrared 
photometry was summarized by Skopal (2000). At positions 
of the inferior conjunction of the giant, the ultraviolet contin- 
uum is Rayleigh scattered by HI atoms of the giant's wind, 
which suggests a high inclination of the orbit (Vogel 1991). To 
show properties of the ultraviolet continuum we selected IUE 
observations taken close to the superior conjunction of the giant 
(SWP42347 + LWP21 103, 28/08/91, <p = 0.47 and SWP55897 
+ LWP31444, 10/09/95, ip = 0.52) and at the opposite position 
(SWP18993 + LWR15045, 13/01/83, <p = 0.95). 

Radiation from the giant. Fluxes of the infrared photometry 
can be compared well with the synthetic spectrum calculated 
for r eff = 3 500K. The scaling factor, k g = 2.23 10~ 17 yields 
the integrated flux F° bs = 1.9 10~ 7 ergcm" 2 s _1 . The timing 
of the eclipse at /11320 A corresponds to R g /A = 0.32 + 0.05 
(see Appendix A of Skopal 2001a). For P or b = 482 days and 
the total mass of the binary, M T = 2.5 Mq (M h = 0.6 Mq 
and elements of Fekel et al. 2000a) we get the separation of 
the stars A = 386%, which yields R g = 124 + 19 Rq. This 
quantity can be converted by 8 g to the distance d = 590 + 92 pc, 
which is well within the range of 400 to 1960pc measured by 
Hipparcos. 

Radiation from the ultraviolet. The last low-resolution 
spectrum of EG And made by IUE (10/09/95) is character- 
ized by a strong nebular emission. Its parameters of our so- 
lution (T e = 25 000 ± 3 000K, k N = 1.25 10 15 cnr 5 ) con- 
strain a minimum temperature of the ionizing source, T™ n 
= 95000K (k b = 1.53 10~ 24 , 8 ~ 1, Eq. 23), which makes 
it just capable of producing the observed emission measure, 
EMobs = 5.2 10 5 V/590pc) 2 cm~ 3 . These parameters im- 
ply L h = 77(c//590pc) 2 L Q , 6 h = 1.24 10~ 12 and R h = 
0. 032(J/ '590 pc) Rq. We adopted this temperature to analyze 
the remaining two spectra. Observations from 28/08/91 show 
a significantly smaller amount of nebular emission (k N = 
2.5 10 14 cm~ 5 ) radiating at a low temperature of T e = 13 500 + 
1 000 K, even though it was taken at similar orbital phase 
as the previous spectrum. However, the luminosity of the 
hot stellar source was comparable: k b = 1.27 10~ 24 (0h = 
1.1 10~ 12 ), which yields R h = 0.029(c//590 pc) R Q and L h = 
60(<f/590pc) 2 Lq. The spectrum taken close to the inferior 
conjunction of the giant (13/01/83) shows a strongly atten- 
uated far-UV continuum by the Rayleigh scattering process 
(«h = 9.4 10 22 cm~ 2 ) and features of the iron curtain absorp- 
tions pronounced at ~ /11400A and ~ /11700A. However, the 
stellar radiation can be scaled with approximately the same 
quantity of k b as on 28/08/91 and 10/09/95, i.e. if one removes 
the effect of the Rayleigh scattering, the hot star luminosity 
remains practically unchanged in all spectra. The nebular com- 
ponent of radiation is fainter at positions near to the giant's 
occultation. Skopal (2001a) showed that the variation in the 
EM follows well that observed in the LCs. Here the ampli- 
tude Am EM = -2.5 log(£ , M m i n /£ , M max ) = 0.36 mag is compara- 
ble with that measured photometrically in the U band (Fig. 2). 
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nebular components of radiation. The solid thick line is the resulting 
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Parameters of the hot star in EG And suggest that the observed 
hot star luminosity is balanced by accretion processes at a rate 
of M acc = a few x lO^Afoyr" 1 (see Skopal 2005, for de- 
tails). Then the observed increase in the nebular emission on 



10/09/95 could be caused by a transient increase in the mass- 
loss rate - an injection of a surplus of emitters into an open 
Hn zone (S < 1 on 28/08/91) increases the rate of ioniza- 
tion/recombination acts. From this point of view, the 0.3 - 
0.5 mag flares in U result from a variable mass-loss rate from 
the giant. 



4.3. Z Andromedae 

Z And is considered as a prototype symbiotic star. The LC is 
characterized by phases of activity with up to 2-3 mag brighten- 
ings, alternating with periods of quiescence (Fig. 3). Recently 
Skopal (2003a) revealed a high inclination of the Z And orbit. 
To show properties of the ultraviolet continuum during quies- 
cence we selected IUE spectra taken at the position when the 
hot star was in front (SWP18601 + LWR14669, 19/11/82, ip 
= 0.49; SWP55160 + LWP30952, 19/11/82, ip = 0.57) and at 
the opposite position (SWP32845 + LWP12624, 03/02/88, <p = 
0.99). For the active phase we selected spectra taken at the 1985 
maximum (SWP27370 + LWP07370, 24/12/85, tp = 0.98). 

Radiation from the giant. We matched the IR fluxes by 
a synthetic spectrum with r e tf = 3 400K and scaling k g = 
2.66 10~ 18 , which corresponds to the bolometric flux, F° hs = 
2.0 10~ 8 ergcm^s -1 . Our value of 9 g is in good agreement 
with that derived from the surface brightness relation for M- 
giants (Dumm & Schild 1998), 6 g = 1.7 10~ 9 , given by the red- 
dening free magnitudes K = 4.85 and J = 6.09 mag. Our value 
of 9 g then gives the radius of the giant R g = 106 (d/ 1 .5 kpc) Rq 
and the luminosity L g ~ 1 400 (d/ 1 .5 kpc) 2 Lq. 

Radiation from the ultraviolet: Quiescent phase. Features 
of the iron curtain can be recognized mainly in the short- 
wavelength part of the spectrum. Our solution for the spec- 
tra at <p ~ 0.5 corresponds to T h = 120 000 K = T™ n 
(see discussion below) and very different electron tempera- 
tures, r e (19/ll/82) = 20500 + 2000K and r e (29/06/95) = 
11500 + 1000K. The significant change in T e by about 
10000K is observationally given by the change of the slope 
of the near-UV continuum profile (Fig. 3). The stellar com- 
ponent of radiation was scaled with k h = 2.73 10~ 24 , which 
yields 6»h = 1.65 10~ 12 , R h = 0.11(J/1.5 kpc) Rq and L h = 
2 300(<f/1.5 kpc) 2 Lq. The nebular component of the radi- 
ation was scaled with %(19/ll/82) = 3.610 15 cirT 5 and 
&n(29/06/95) = 2.7 10 15 cm -5 . These fitting parameters give 

5 = 0.92 and 1.1 on 19/11/82 and 29/06/95, respectively. The 
case of 6 ~ 1 corresponds to a lower limit of the hot star 
temperature, T™" ~ 120 000 K, under which the hot star is 
not able to produce the observed amount of nebular emission 
(Eq. 23). However, parameters of Z And suggest a very open 
Hn zone (e.g. Fig. 7 of Fernandez-Castro et al. 1988), which 
means that a certain fraction of the ionizing photons escapes 
the system without being converted to nebular radiation (i.e. 

6 <sc 1). This situation requires even higher at the same 
scaling to balance the measured EM. Therefore the real val- 
ues of Lh and R^ are larger/smaller than we derived from our 



model for the T™ 1 " temperature. The spectrum from 03/02/88 
(<p ~ 0) is characterized by the Rayleigh attenuated far-UV con- 
tinuum with n H ~ 3.9 10 22 cirT 2 and features of the iron curtain 
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Fig. 3. As in Fig. 2, but for Z And during quiescent and active phases. 
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that are more pronounced than at the opposite binary position. 
This can be interpreted in terms of the atmospheric eclipse due 
to Rayleigh scattering of the far-UV photons on the hydrogen 
atoms of the neutral giant's wind (Isliker et al. 1989). During 
quiescent phases this effect is a strong signature of a highly in- 
clined orbit of the binary. A further distinct difference between 
the spectra taken at <p ~ 0.5 and ip ~ is an attenuation of 
the UV continuum at all wavelengths by a factor of about 2. 
This effect is present in all IUE spectra of Z And in the range of 
phases (p « + 0. 15. As a result such an extinction process re- 
sults in lower values of the derived parameters, 9^ < 1.0 10~ 12 , 
Rf < 0.068(,i/ 1 -5 kpc) R Q and L h > 860(c// 1 .5 kpc) 2 Lq. Also 
the amount of the nebular emission decreased by a factor of 
about 5 and the parameter 6 = 0.2 for = 1.0 10~ 24 , which 
reflects a significant attenuation/loss of the nebular photons in 
the direction of the observer at ip ~ 0. 

Radiation from the ultraviolet: Active phase. The selected 
spectrum is characterized by a strong Rayleigh attenuation of 
the far-UV continuum and the effect of the iron curtain absorp- 
tion, which drastically depressed the level of the continuum in 
between 2150 and 3 000 A (Fig. 3). Note that during activity 
these effects are also significant out of the occultation by the 
giant (e.g. on SWP22684 + LWP03099, <p = 0.15). The tem- 
perature of the stellar component of the radiation significantly 
decreased to ~ 19 000K and is transferred throughout the 
~ 2 10 23 crcT 2 hydrogen atoms on the line of sight. Below we 
call this stellar component of radiation during active phases the 
'hot stellar source' (HSS). There are two components of neb- 
ular radiation: (i) A low-temperature nebula (LTN) with typi- 
cal properties of quiescent phases and (ii) a high-temperature 
nebula (HTN) corresponding to a high electron temperature 
(T e ~ 40 000 K). The latter is required to fit the far-UV contin- 
uum. With analogy to AR Pav (see. Skopal 2003b), we scaled 
a 40 000 K hot nebular radiation of hydrogen to the non-zero 
continuum around the Lya line. Ignoring the effect of the iron 
curtain leads to an artificially high He ++ abundance and too low 
electron temperature of the LTN (see Fig. 3 of Skopal 2003a). 



4.4. AE Ara 

Visual LC of AE Ara (1987- 2004) displays some bright stages 
alternating with periods of quiescence (Fig. 4). For the purpose 
of this paper we used the only well exposed low-resolution 
spectrum from the IUE archive (SWP06793+LWR05784, 
08/10/79, <p = 0.53). To flatten its pronounced 2 200 A feature 
we determined the colour excess, £ B -v = 0.25 +0.05. 

Radiation from the giant. A synthetic spectrum with T e g = 
3 200K matches best the IR fluxes from IJKL photometry 
(Allen & Glass 1974; Munari et al. 1992). Its scaling, 
k g = 8.1 10~ 19 , yields F° hs = 4.8 10~ 9 erg cm" 2 s" 1 . If we 
adopt the radius of the giant, R g = 130 - 150 Rq, in agree- 
ment with the empirical relations between the spectral type 
and T e g (e.g. Belle et al. 1999), then 6 g is transformed to 
d = 3.5 + 0.3 kpc, where the uncertainty is given only by the 
adopted range of the giant's radius. This result is in excellent 
agreement with that derived from the surface brightness rela- 
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Fig. 5. The SED for CD 43° 14304. It is very similar to that from active 
phases of other well known systems. Crosses (x) represent fluxes in 
the continuum of the optical spectroscopy carried out in July 1991 
(Gutierrez-Moreno et al. 1999). At the time of IUE observations the 
system just entered an active phase (top panel). 



tion for M-giants (Dumm & Schild 1998). The giant's lumi- 
nosity L g = 1 800 + 300 (J/3.5 kpc) 2 L (Eq. 4). 

Radiation from the ultraviolet. The hot star in AE Ara dur- 
ing the time of the IUE observation was in an active stage 
with a significant contribution in the optical. A moderate in- 
fluence of the iron curtain absorptions can be also recognized. 
Our model identified a warm pseudophotosphere (the HSS) 
radiating at T h = 23 000 + 10000/ - 3 000K. The scaling 
factor, &h = 1.40 10~ 22 , implies its effective radius, = 
1.8(c//3.5 kpc)/v0 (Eq. 6). To obtain a better fit of the long- 
wavelength end of the spectrum, a nebular radiation determined 
byr e = 17 000 + 3 000 K and the scaling, k N = 3.6 10 14 ctrT 5 , 
was required. There is no direct signature for the presence 
of the HTN component in the spectrum. The active stage of 
AE Ara is signaled by its two-temperature UV spectrum, which 
develops during active phases of all the objects with a high or- 
bital inclination (Sect. 5.3.4). 



4.5. CD-43 14304 

Visual estimates of this system were regularly carried out by 
Albert Jones (Fig. 5, private communication). The LC shows 
periods of a lower and higher level of activity differing by 
about 1 mag. Only a few multicolour photometric observations 
were made by Munari et al. (1992) and Gutierrez-Moreno 
et al. (1999) in the optical UBVRI and infrared JHK pass- 
bands. Schmid & Nussbaumer (1993) found that the light of 
CD-43° 14304 is reddened with £ B -v < 0.2. The observed in- 
dex V - K = 4.0 and that corresponding to the 4 000 K giant, 
(V - K) g = 3.6 (Belle et al. 1999), yields £ B -v = 0. 15. This 
value represents a lower limit, because of a contribution from 
other sources in the V band. Therefore we adopted £b-v = 0.2. 
There is only one spectrum in the final IUE archive (SWP41632 
+ LWP20373, 15/05/91, <p = 0.95) applicable to our modeling. 
The spectrum taken by the long wavelength prime was poorly 
exposed (Fig. 5). 
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Radiation from the giant. The (B)VRIJK photometric mea- 
surements can be compared well with a synthetic spectrum 
of r eff = 4000K scaled with k g = 1.9 10 19 , which gives 
F° hs = 2.7 1(T 9 ergctrT 2 s _1 . If we adopt the giant's radius to 
be R g = 40 Rq as a typical value for the spectral type K7 III 
(Murset & Schmid 1999; Belle et al. 1999) then the angular 
radius yields d = 2.1 kpc. Accordingly, the giant's luminosity 
L g = 380(d//2.1kpc) 2 L Q . 

Radiation from the ultraviolet The two-temperature type 
of the ultraviolet spectrum signals an active stage of CD- 
43°14304. Independently, Gutierrez-Moreno et al. (1999) 
found an increase of the intensity shortward of 5 000 A in their 
optical spectra from 1991 June and July (i.e. taken near to the 
IUE observations) with respect to the spectrum from 1987 by 
a factor of 3 at 4 000 A. They ascribed this optical brightening 
to a mild outburst. Our solution indicates a cool shell (i.e. the 
HSS) re-radiating the hot star photons in the direction of the 
observer at T h = 20000 + 7 000/ - 3 000K and scaling of 
feh = 2.4 10 -23 , which results in a significant contribution in the 
optical. The spectrum shows the far-UV continuum Rayleigh 
attenuated by 2.5 10 23 cuT 2 H i atoms and pronounced features 
of the iron curtain absorption. The position of the binary at the 
time of observation (<f = 0.95) enhances both effects. The ef- 
fective radius of the shell is R 1 ^ = 0.46(<f/2.1 kpc)/?0 and its 
luminosity Lh = 30(c//2.1 kpc) 2 L© CEq. 6). A non-zero level 
of the Rayleigh attenuated far-UV continuum at/around Lya in- 
dicates a contribution from the HTN with k N = 1.9 10 14 cirT 5 
for T e = 40 000 K. The continuum profile at the end of the 
LWP spectrum and the following part containing the Balmer 
jump reflects the presence of a faint radiation from the LTN: 
T e ~ 12000 K, k N ~ 0.5 10 14 cirT 5 . Evolution in the Ha line 
profile along the orbital motion (Schmid et al. 1998) is very 
similar to that observed in eclipsing systems. For example, 
EG And (Oliversen et al. 1985), V1329Cyg (Ikeda & Tamura 
2000), AXPer(Iijima 1988), RW Hya (Schild et al. 1996) and 
SYMus (Schmutz et al. 1994). This suggests that the orbital 
plane of CD-43° 14304 is highly inclined to the observer, which 
is consistent with the presence of the two-temperature ultravi- 
olet spectrum. 



1900 



1920 



1940 



1960 



1980 



2000 



10 



TX CVn 

B a 

V + 



L_. .. ; \ y i t t l ^» 



IUE 



20000 



30000 40000 
Julian date - 2 400 000 



50000 





2.5 - 




TX CVn 


T eff =4100K 


+ 


2 - 
1.5 - 




Active phase 


T h = 17000 K 

t 6 = ioooo k ; 


'< 


SWP15655+ 


05/12/1981 


n H = 2.3E+22 cm" 2 




1 - 


+LWR12083 






u 


0.5 - 








,03, 





i ip 






g Flu? 


-0.5 - 




'jlljl 




o 


-1 - 




U B V H I 


. J H K L Tk 






3.2 3.4 


3.6 3.8 4 


4.2 4.4 4.6 



TX CVn 




SWP05677: 29/06/79: <p=0.57 - - 
SWP15655+LWR12083: 05/12/81 : <p=0. 06 



t~nr. 




3.3 3.4 3.5 

log Wavelength [A] 
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4.6. TX Canum Venaticorum 

The historical LC shows the eruptive character of TX CVn. The 
present active phase is maintained from the last major outburst 
in 1962 (Fig. 6). Recently, Skopal et al. (2004) revealed two 
minima in the U-LC, whose timing coincides with that of the 
inferior conjunction of the giant (Kenyon & Garcia 1989). 
This suggests a high orbital inclination of TX CVn. Here we 
used well-exposed IUE spectra taken close to the inferior and 
superior occultation of the giant (SWP15655 + LWR12083, 
05/12/81, <p = 0.06 and SWP05677, 29/06/79, <p = 0.57). 

Radiation from the giant. A synthetic spectrum of = 
4100K matches well the VRJHKL photometric measure- 
ments. Its scaling, k g = 4.5 10~ 19 , corresponds to F° bs = 
7.1 10~ 9 erg cirT 2 s _1 . According to the empirical relations be- 
tween the spectral type, V - K index and (e.g. Belle et al. 
1999), the giant's radius R g = 30 + 10 Rq, which then trans- 



forms the 9 g parameter to the distance d = 1.0 + 0.3 kpc. The 
corresponding luminosity L g = 230 (d/ 1.0 kpc) 2 Lq. 

Radiation from the ultraviolet. The profile of the ultraviolet 
continuum displays characteristics of an active phase in agree- 
ment with evolution in the LC (Fig. 6). The effect of Rayleigh 
scattering and absorptions of the veiling gas are present at very 
different orbital phases (</> = 0.06 and 0.57), which suggests its 
distribution around the central star. However it is not possible to 
determine its geometry more accurately, because the spectrum 
does not contain emission lines - the shell probably encom- 
passes the whole star core (see also the discussion on the nature 
of the hot component by Kenyon & Garcia 1989). Our solu- 
tion for the 05/12/81 spectrum suggests a black-body radiation 
at T h ~ 17 000K (k h = 2.6 10~ 22 ) as a dominant component 
from the hot object. Only a very small nebular contribution is 
required to fit better the long-wavelength end of the spectrum. 
However, the line blanketing affects drastically the ultraviolet 
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Fig. 7. The UV/optical/IR SED for TCrB. A black-body accretion disk 
dominates the UV domain with a significant contribution in UB pass- 
bands. The Rayleigh scattering and absorptions of the iron curtain 
affect the UV continuum. The contribution from the nebula is very 
small. Dates of observations are marked in the LC at the top panel. 



continuum, which does not allow us to determine reliably the 
continuum fluxes and thus to distinguish unambiguously the 
presence of a faint LTN in the spectrum. 

4.7. TCoronae Borealis 

TCrB is a symbiotic recurrent nova with two major outbursts, 
recorded in 1866 and 1946. Recent evolution in the LC shows 
periods of higher and lower states of activity (e.g. Stanishev 
et al. 2004, Fig. 7). Currently, it is believed that the hot star 
luminosity is powered by an accretion process triggered by a 



mass transfer from the lobe-filling giant. Selvelli et al. (1992) 
suggested a white dwarf as the accretor to satisfy mainly the ul- 
traviolet IUE observations, while Cannizzo & Kenyon (1992) 
tested accretion of a gaseous torus onto a main-sequence star to 
explain the secondary outburst of the 1946 eruption. To model 
the ultraviolet continuum we selected observations from the gi- 
ant's inferior conjunction (SWP31095 + LWP10902, 03/06/87, 
^ = 0.99) and from the opposite site (SWP19869, 01/05/83, 
ip = 0.42 and LWR16145, 13/06/83, tp = 0.62 - no simultaneous 
SWP/LWR exposures are available at this position). 

Radiation from the giant The red giant's fluxes can be 
matched by a synthetic spectrum of r eff = 3 400K and the 
scaling, k g = 3.13 10~ 18 , which corresponds to the observed 
bolometric flux, F° bs = 2.37 10~ 8 ergcm^s -1 and the angu- 
lar radius, 6 g = 1.77 10~ 9 . These parameters imply the gi- 
ant's radius, R g = 75 + 12(<f/960pc)/?o and the luminosity, 
L g = 680 + 230(c//960pc) 2 L Q (Table 1). 

Radiation from the ultraviolet. Both spectra display the 
Rayleigh attenuated far-UV continuum and pronounced fea- 
tures of the iron curtain. This constrains the veiling material 
surrounding the hot source, being located at/behind the outer 
rim of the disk, extended vertically from the disk plane, so as 
to cause the observed absorptions. At <p ~ the effect of the ab- 
sorbing slab is enhanced by the cool giant's wind. We estimated 
flux-points of the continuum with the aid of calculated spectra 
including the blanketing effect (Shore & Aufdenberg 1993). 
We performed modeling of the continuum by substituting the 
single Planck function in Eq. (13) by a component of radiation 
from an accretion disk. Then Eq. (13) in the ultraviolet is 

= k AD x F A (AD) e-"^ +k N x £i(H, He + , T e , a). (24) 

The function F ^(AD) represents a flux distribution from an op- 
tically thick accretion disk that radiates locally like a black 
body. Here we adopted expressions according to the textbook 
of Warner (1995), recently used also by Skopal (2003b). A 
fundamental parameter in models with an accretion disk is the 
disk temperature r* = 2 r max , which determines the slope 
of the UV continuum. Observations from 03/06/87 suggest 
r t * 1 10 000 K at all wavelengths. We found that it is possible 
to match the selected points of the continuum by very different 
sets of model parameters. Below we summarize all possible 
solutions to select ranges of fitting parameters that are physi- 
cally plausible. In all calculations we adopted the core mass of 
the accretor to be 1 .37 Mq and the orbital inclination i = 67° 
(Stanishev et al. 2004). 

1 . Solutions with the accretion disk only: 

(a) If the radius R^ ~ /? W d (we adopted R^ = 0.01 Rq) then 
the temperature required by the slope of the continuum, ~ 
1 10 000 K, corresponds to the mass accretion rate M acc ~ a few 
xlO~ 9 M0yr~', which yields L acc ~ a few xL@. These values 
are much lower than those given by observations, L \, S (AD) = 
J A F A (AD) dA I cos(;) ~ 300 Lq. 

(b) If the radius R^ > /?wd (we tested a range of 0.10 
to 0.05 Rq), the temperature can be satisfied by M acc ~ 
10~ 6 Mq yr~' corresponding to L acc ~ L Q b s (AD). 

2. Solutions with the accretion disk and nebula. These mod- 
els are not unambiguous. Models with a higher express bet- 
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ter only the far-UV region, which then requires a strong nebular 
contribution to fit the near-UV part of the spectrum. 

(a) The case of ~ R WD and a large nebular contribu- 
tion leads to » 110 000 K, L acc ~ 10 3 L© and M acc ~ 
1-10 10~ 7 Mq yr -1 . These models are not consistent with ob- 
servations in too high L acc and the nebular emission. Models 
with small contributions from the nebula converge to the cases 
of point 1 (a) above. 

(b) Solutions for Rh > Rwd correspond to ~ 1.1 — 
1.3 10 5 K, M acc ~ 1 - 1.6 X 10- 6 M o yr 1 and L acc ~ L obs (AD) ~ 
300 - 400 Lq. The emission measure is small and it is not pos- 
sible to determine T e unambiguously. 

Our modeling of the SED of T CrB shows that a major frac- 
tion of the UV radiation is produced by an optically thick ac- 
cretion disk. Models whose parameters are consistent with ob- 
servations result from an accretion process at high accretion 
rates of M acc ~ 1 - 1.6 x 10~ 6 Mq yr~' onto the accretor with 
Rh > Rwd (Fig. 7). These results, however, do not support the 
presence of the white dwarf as the accretor in TCrB. 

4.8. BFCygni 

BFCyg is an eclipsing symbiotic star (Skopal 1992, Fig. 8 
here). Its historical LC shows three basic types of active 
phases - nova-like and ZAnd type of outbursts and short- 
term flares (Fig. 1 of Skopal et al. 1997). Here we selec- 
ted observations from the optical maximum (SWP39163 + 
LWP18251, 30/06/90, cp = 0.59), the eclipse effect (SWP41531 
+ LWP20275, 01/05/91, ip = 0.0), a transition to quiescence 
(SWP44370 + LWP22778, 10/04/92, ip = 0.46) and a quies- 
cence (SWP58385 + LWP32695, 29/09/96, <p = 0.61). 

Radiation from the giant. Fluxes of the photometric mea- 
surements can be compared with a synthetic spectrum of r e ff = 
3 400K (Figs. 8). Its scaling, k g = 4.5 10~ 19 , corresponds to 
the bolometric flux, F° bs = 5.9 10~ 9 ergcirT 2 s~'. Empirical 
relations between the spectral type and T e g suggest the giant's 
radius R g ~ 150 Rq, which yields the distance d ~ 3.8 kpc. 
On the other hand the contact times of the 1991 eclipse give 
Rg/A = 0.54 + 0.02, which suggests a very large radius of 
the giant (Skopal et al. 1997). A solution for the spectro- 
scopic orbit (Fekel et al. 2001) yields the Roche lobe ra- 
dius Rl/A = 0.49. These results imply the lobe-filling giant 
in BFCyg. However, the broad eclipse profile may be partly 
caused by attenuation of the optical light in the extended dense 
atmosphere of the giant. Note that the SWP42356 spectrum 
from 29/08/91 shows a steep far-UV profile, which signals 
the out-of-eclipse position, while the brightness in the opti- 
cal remains at the level of the total eclipse. Therefore here we 
adopted R g = 150 Rq. 

Active phase. The continuum is drastically depressed by the 
iron curtain absorptions and the far-UV continuum profile in- 
dicates a hydrogen column density with «h ~ 5 10 22 crrT 2 and a 
contribution from the HTN. Other spectra from the maximum 
are of the same type, i.e. its profile does not depend on the or- 
bital phase. Only very faint emissions of Siiv 1403, Hen 1640 
and Oiii 1664 A can be still recognized. Our solution suggests 
the presence of a strong stellar component (r n ~ 21 500 K, ~ 
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1.8 10~ 21 ) superposed on the HTN emission (T e ~ 42 000 K, 
k N ~ 6.9 10 15 cirT 5 - as derived from the following eclipse). 
Both components dominate the optical and in part influence 
also the IR wavelengths. This composition explains colour in- 
dices during the totality of the 1991 -eclipse, which also reflect 
a combined spectrum (Skopal 1992). 

Eclipse. The ultraviolet continuum from the 1991 -eclipse is 
flat in profile radiating at a level of ~ 10~ 12 ergcrrT 2 s _1 A -1 . 
We ascribe this component to the emission from the HTN 
(r e ~ 42000K, EM mN ~ 1.1 10 61 (W/3.8kpc) 2 cm- 3 and 
Lhtn ~ 2 100 (J/3.8 kpc) 2 L©). 

Transition to quiescence. First we subtracted the HTN con- 
tribution and fitted the rest fluxes adopting the hot star temper- 
ature T h = 95 000 K = r™ n (6=1, Eq. 23). The presence of 
the HTN emission is given by observations from the preceding 
eclipse. A high 7h here is required by a large emission measure 
EM = 1.210 61 (J/3.8kpc) 2 cm 3 in total (Table 3). We note 
that a fraction of the nebular emission can be produced by col- 
lisions as suggested by the presence of the HTN (Sect. 5.3.1). 
This would allow the ionizing source to be cooler than T™ n . 
Iron curtain features are relatively faint and the far-UV spec- 
trum is attenuated with only n H ^ 10 22 cm~ 2 . If one ignores in- 
fluence of the iron curtain a formal good fit with 7h = 60 000 K 
is also possible (Fernandez-Castro et al. 1990). 

Quiescent phase. The stellar component has practically 
the same properties as that from the transition period, but 
without signatures of more pronounced absorptions from the 
veiling neutral gas. However, the nebular component (r e = 
22 800 K) decreased by a factor of about 4, to EM = 
3.1 10 60 (<f /3.8 kpc) 2 cm -3 , despite of that both the spectra were 
taken at the same orbital phase (tp « 0.5). A decrease of the 
emission measure at the same ionizing capacity of the hot star 
implies that a fraction of the ionizing photons escapes the sys- 
tem without being converted into the nebular emission, i.e. the 
parameter 6 decreased and the nebula became more open than 
during the transition period. A constant mass-loss rate from the 
giant would mean that the hot star wind, which developed dur- 
ing the activity, caused a surplus of ionizations. In part it is 
itself a subject of ionization and in part it creates a fraction of 
the HTN by collisions. Otherwise the mass-loss rate from the 
giant had to significantly decreased. 

4.9. CHCygni 

CHCyg is the brightest symbiotic star in the infrared pass- 
bands (J ~ 0.9, K ~ -0.6 mag). It is at 270 + 66 pc 
(Viotti et al. 1997). We consider CHCyg to be a triple- 
star system consisting of the inner, 756-day period symbi- 
otic binary, which is moving in a common 5300-day period 
outer orbit with another red giant (Hinkle et al. 1993; Skopal 
et al. 1996). To demonstrate evolution in the SED during 
active phases of CHCyg we selected observations from its 
1992-95 activity: Before the 1994 eclipse in the inner binary 
(SWP47771 + LWP25627, 30/05/93, (p = 0.30), during the to- 
tality (SWP52989 + LWP29646, 05/12/94, <p = 0.03) and after 
the eclipse (SWP54254 + LWP30341, 30/03/95, <p = 0.18). 
Radiation from the giant. The slope of the VRIJ fluxes 
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Fig. 9. Example of the SED for CHCyg during its 1992-95 active 
phase. Spectra prior to, in and after the 1994 eclipse were selected. 
Corresponding dates are marked at the LC on the top panel. 



suggests a very low temperature of the giant(s) in the system. 
We matched the U BV fluxes from the giant and the observed 
fluxes in the RIJKLM bands by a synthetic spectrum of r e tf = 
2600K (see Appendix B of Skopal et al. 2002a, for details). 
Corresponding parameters, F° bs = 2.4 10 _6 ergcm _2 s _1 and 
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g = 3.1 1(T 8 yield the giant's radius R g ~ 370(<f/270pc)/? o 
and the luminosity L g ~ 5 600 (<i/270pc) 2 Lq. These parame- 
ters probably belong in major part to the giant at the outer orbit 
in the triple-star model. We note that the presence of the two 
red giants in the system has been suggested by the timing of 
photometric eclipses (e.g. Skopal 1997) and by the COAST 
interferometric measurements at 0.9 pm (Young et al. 2000). 

Eclipse. The spectrum taken during the totality is almost flat 
in profile with a moderate increase of fluxes from about 2 000 A 
to the long-wavelength end of the spectrum (Fig. 9). On the 
radio maps, Crocker et al. (2001) found two components of 
nebular emission. One is of thermal nature located around the 
central bright peak and another one is of non-thermal nature 
associated with the extended regions. Eyres et al. (2002) sug- 
gested that the extended emission has electron temperatures of 
~ 100 000 K. Therefore we compared a 10 5 K hot nebular radia- 
tion to the SWP spectrum and filled in the rest of the IUE spec- 
trum by a 10 4 K nebular radiation. The resulting fit expresses 
satisfactorily the overall profile. The long -wavelength part of 
the spectrum is affected by the iron curtain absorptions, which 
is consistent with the location of the LTN at the central bright 
peak on the HST image. 

The out-of-eclipse SED. During active phases the profile of 
the out-of-eclipse continuum is of the same type (here Fig. 9 
and Figs. 1, 2 of Skopal et al. 1998). A relatively cool pseu- 
dophotosphere (here 7h ~ 14 000 K) is drastically affected by 
absorption of the veiling gas and the far-UV continuum is at- 
tenuated by a few xl0 22 cm~ 2 hydrogen atoms. These effects 
are present at all orbital phases. Modeling the continuum of the 
30/05/93 spectrum we neglected the nebular component of ra- 
diation due to a very low level of the continuum at Lya. Finally, 
we note that including the influence of the iron curtain in mod- 
eling the ultraviolet continuum changes the interpretation sug- 
gested by Skopal et al. (1998). 

Quiescent phases. All components of radiation associated 
with the hot star are negligible during quiescent phases (cf. 
Fig. 1 of Skopal et al. 1998). Such a behaviour represents a 
significant difference from all other symbiotic stars with rele- 
vant observations. 



4.10. CICygni 

The last active phase of CICyg began in 1975 (Belyakina 
1992). During the first four cycles from the maximum, nar- 
row minima - eclipses - developed in the LC; a typical fea- 
ture of active phases of symbiotic stars with a high orbital 
inclination. From about 1985 the profile of minima became 
very broad indicating a quiescent phase (Fig. 10). Accordingly 
we selected observations from the active phase (SWP03816 + 
LWR03396, 05/01/79, (p = 0.39), from the following eclipse 
(SWP09255 + LWR08003, 10/06/80, ip = 0.0) and the quies- 
cence (SWP36321 + LWP15571, 22/05/89, <p = 0.82). 

Radiation from the giant. Fluxes corresponding to the IJKL 
photometry can be matched by a synthetic spectrum of T e g = 
3 300K scaled with k g = 4.0 10~ 18 (6 g = 2.0 10~ 9 ), which 
yields the integrated flux, F° hs = 2.8 10~ 8 ergcm _2 s _1 . We 
note that the surface brightness relation for M-giants gives a 
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Fig. 10. Examples of the SED during the activity (A), eclipse (E) and 
quiescence (Q) for CI Cyg. Corresponding dates are marked in the LC. 



similar value of 9 g =2.1 10~ 9 for the observed reddening-free 
magnitudes, K = 4.37 and J = 5.63 mag. Timing of the eclipse 
profile, which yields R g /A = 0.38+0.02, and other fundamental 
parameters derived by Kenyon et al. (1991) suggest the giant's 
radius R g ~ 180 Rq. This implies the distance d ~ 2.0 kpc and 
the luminosity L g ~ 3 400 (c//2.0kpc) 2 Lq. 
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Radiation from the ultraviolet: Active phase. The profile of 
the UV continuum during the active phase of CI Cyg is of the 
same type as that of other active symbiotics with a high or- 
bital inclination - Rayleigh attenuated radiation from a warm 
HSS, which is affected by the iron curtain absorptions (here 
T h ~ 28 000 K, k h ~ 5.7 10~ 23 and n H ~ 1 10 23 crrT 2 ). This 
component represents the main contribution in the UV/U spec- 
tral region. In addition, two components of the nebular radi- 
ation, from the HTN and LTN, as described for ZAnd, are 
required to get a satisfactory fit. The HTN emission is seen 
directly during the eclipse as the only component and out of 
the eclipse its presence is signaled by the non-zero level of the 
Rayleigh attenuated far-UV continuum (Fig. 10). 

Radiation from the ultraviolet: Quiescent phase. In fit- 
ting the ultraviolet continuum we fixed the hot star tempera- 
ture to its Zanstra temperature, T h = 115 000 K (MNSV, see 
Sect. 5.2.2 for details). Our solution yields a rather high elec- 
tron temperature for the nebula, T e = 24 000 ± 4 000 K, which 
is constrained by a flat profile and a high level of the long- 
wavelength part of the continuum (A > 2 000 A). The scaling 
factor of the nebular component, fc N = 1.0 10 15 crrT 5 , gives 
EM ohs = 4.8 10 5 V/2.0kpc) 2 cirr 3 . 



4.11. V1329Cygni 

VI 329 Cyg is a symbiotic nova which erupted in 1964 and 
peaked at m pg = 11.5 mag in 1966 October. Prior to this 
outbreak, VI 329 Cyg was an inactive star of about 15th mag- 
nitude displaying ~2mag deep eclipses (see Fig. 1 of Munari 
et al. 1988). After the outburst the symbiotic phenomenon 
with characteristics of a classical symbiotic star developed. 
For the purpose of this paper we selected observations from 
the maximum (SWP19577 + LWR15609, 29/03/83, (p = 0.57) 
and minimum (SWP14688 + LWR11266, 07/08/81, cp = 0.95) 
of the wave-like variation in the optical LCs (Fig. 11). To 
better visualize their profile we smoothed them by averaging 
the measured points within 10 A bins. 

Radiation from the giant. Infrared fluxes were derived from 
average values of RJKL measurements published by Tamura 
(1983), Lorenzetti et al. (1985), Taranova & Yudin (1986), 
Munari et al. (1988), Nussbaumer & Vogel (1991) and Munari 
et al. (1992). They can be matched by a synthetic spectrum 
of r eff = 3 300K scaled to F° hs = 3.4 10~ 9 erg cnr 2 s -1 , 
(k s = 5.2 10~ 19 , 6 g = 7.2 10- 10 ). The same value of 6 g can 
be obtained for the M-giant's magnitudes, K = 6.67 and J 
= 7.88 mag, by using the Dumm & Schild (1998) relation. 
Combining these parameters with those recently derived 
by Fekel et al. (2001) (Af g = 2.2 Mq, M h = 0.75 Mq, 
R g = 132 ± 40 Rq) we get the distance d = 4.2 ± 1.3 kpc, and 
the luminosity L g = 1 900 (c//4.2kpc) 2 Lq, where the uncer- 
tainty reflects only that in R g . On the other hand, the maximum 
width of the pre-outburst eclipses, U — t\ ~ 0.13 P or b (Fig. 1 
of Schild & Schmid 1997), yields R g ~ 228% assuming 
the eclipsed object to be a point source. However, such a large 
radius of the giant could result from an atmospheric extinction 
as in the case of BFCyg (Sect. 4.8). Therefore we prefer the 
value of 132 Rq, which was derived under the assumption that 
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Fig. 11. The SED in the continuum of VI 329 Cyg at maximum and 
minimum of the wave-like variation in the LC. Crosses (x) at the SED 
from the maximum represent continuum fluxes taken from the optical 
spectra presented in Figs. 192 and 193 of Munari & Zwitter (2002). 
Dates of the IUE observations are marked in the LC. 



the giant rotates synchronously with the orbital revolution. 

Quiescent phase. Radiation of V1329Cyg in the 
UV/optical domain is characterized by a very strong neb- 
ular emission (T e = 18 000 K, Ic N = 1.1 10 15 ctrT 5 ), which 
dominates the spectrum in the range of 1 660 - 7 000 A at the 
maximum (Fig. 1 1). This constrains the very high temperature 
of the ionizing source, T h > 170000K = T™" scaled to the 
far-UV continuum with kh = 2.91 10~ 25 , to produce the ob- 
served emission measure, EM b s = 2.3 10 60 (<f/4.2kpc) 2 cirr 3 
at the maximum (i.e. to get the parameter 6 < 1, Eq. 23). 
We note that T™ 1 " required by our model is higher then the 
Zanstra temperature of 145 000 K (MNSV). At the minimum 
the far-UV continuum faded by a factor of about 2 as in the 
case of ZAnd (Sect. 4.3) and the nebular emission decreased 
by a factor of nearly 5 (Table 3). The latter is responsible for 
the large amplitude of the periodic variation in the LC. 
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Fig. 12. Approximate SED for LT Del (LWP spectrum was not avail- 
able). The SWP spectrum was taken at the maximum of the wave- 
like brightness variation in U. Full circles in the SED (•) are photo- 
metric fluxes from the maximum, while asterisk (*) from the mini- 
mum. Crosses (x) represent continuum fluxes from the optical spec- 
trum (Munari & Buson 1992). Small values of the optical/IR fluxes 
from the giant suggest that LTDel is a rather distant object. 



4.12. LTDelphini 

This yellow symbiotic star shows very pronounced orbitally- 
related brightness variation in the U band. This suggests a 
strong and variable nebular emission in the system (Fig. 12). 
In 1994, Archipova et al. (1995) detected for the first time an 
optical outburst of this star. There is only one spectrum in the 
IUE archive taken by the short-wavelength prime (SWP39089, 
14/06/90, (p = 0.49). In spite of this we included this object 
in our sample to get a better representation for S-type yellow 
symbiotics. We determined a new value of £b-v = 0.20 by red- 
dening 1 - 2 10 5 K Plank functions to fit the observed contin- 
uum points of the SWP spectrum. Note that the amount of the 
interstellar extinction at the far-UV is similar to that at 2 200 A. 
A previous value of £b-v = 0.35 (Munari & Buson 1992) led 
to a strong enhancement of the far-UV fluxes, far beyond of 
any reasonable model. 

Radiation from the giant. Fluxes of VRIJKL photomet- 
ric measurements can be matched by a synthetic spectrum of 
r eff = 4100K and scaling of k g = 3.0 lO" 20 , which yields 



the bolometric flux, F° bs 



4.8 10 10 erg cm 2 s Our value 



of r e ff is markedly lower than that suggested by the current 
spectral classification as K0 (Miirset & Schmid 1999). With 
respect to this disagreement we note that spectral classifica- 
tion of yellow symbiotics is based on features of G band, 
Ca 1/14227, Fe l14405 and absorptions of Fei, Srn and the CH 



band. Miirset & Schmid (1999) pointed out that this approach 
may be inaccurate for stars with peculiar abundance patterns 
(here AG Dra, LT Del, BD-21°3873). Our effective temperature 
from IR fluxes suggests a later spectral type, K3 - K4, of a nor- 
mal giant (Belle et al. 1999). Consequently we adopted the 
corresponding average giant's radius of R g ~ 30 Rq, which 
gives the distance d = R g /0 g ~ 3.9 kpc and the luminosity 
L g = 230 (<f /3. 9 kpc) 2 Lq. We note that the giant's radius of 
30% is well inside the range of 18 to 36(M g /1.5 Mq) 1/2 R q , 
which can be derived from atmospheric log(g) = 1.8 + 0.3 
(Pereira et al. 1998). 

Radiation from the ultraviolet. Having only SWP spec- 
trum and a few flux-points from the simultaneous optical spec- 
troscopy (Munari & Buson 1992) the resulting model is ap- 
proximate, mainly in T e with an uncertainty of +3 000 K. = 
100 000 K was adopted as a characteristic temperature. The 
model is shown in Fig. 12. 



4.13. AGDraconis 

The LC of AG Dra shows numerous brightenings by 2-3 mag 
in U with maxima often separated approximately by 1 year. 
Recently, two major eruptions were recorded during 1981-83 
and 1994-96 accompanied by many short-term events lasting 
from a few weeks to a few months. Gonzalez-Riestra et al. 
(1999) identified cool (1981-83 and 1994-96) and hot (e.g. 
1985-86) outbursts differing in their Zanstra temperatures. The 
quiescent LC is characterized by a periodic wave-like varia- 
tion. The top panel of Fig. 13 demonstrates these character- 
istics. The data are from Skopal et al. (2004) and Leedjarv 
et al. (2004). There are no signs either in the optical or far- 
UV regions of eclipses. Schmid & Schild (1997) derived the 
orbital inclination ;' = 60(+8°.2). To demonstrate the nature 
of the observed light variation we selected two IUE spectra 
from quiescence taken at different spectroscopic conjunctions 
of the binary components (SWP37473 + LWP16675, 27/10/89, 
y = 0.62; SWP06650 + LWR05691, 25/09/79, <p = 0.91) and 
two spectra from active phases; one from the major 1994 cool 
outburst (SWP51632 + LWP28752, 28/07/94) and one from 
the hot 1985 eruption (SWP25443 + LWP05513, 13/03/85). 

Radiation from the giant. VRI measurements from the light 
minima during quiescence (recently Leedjarv et al. 2004) and 
the infrared JHKLM photometry of AG Dra can be matched by 
synthetic spectra with T e $ = 4 100 - 4 300 K. According to the 
analysis of Smith et al. (1996) we adopted here T eS = 4 300 K. 
The scaling factor of this spectrum, lc g = 4.9 10~ 19 , corre- 
sponds to the bolometric flux F° hs = 9.51 10~ 9 ergcm _2 s _1 . 
From the surface gravity, log(g) = 1.6 + 0.3 (Smith et al. 
1996), one can derive the giant's radius as R g = (33 + 
ll)(M g /1.5M ) 1/2 fl o and the distance d = R g /6 g = (1.1 + 
0.4)(M g /1.5MQ) 1/2 kpc, where the giant's mass of 1.5 Mq was 
taken from Mikolajewska et al. (1995). The giant's luminosity 
L g = 360 (d/ 1.1 kpc) 2 L Q . 

Radiation from the ultraviolet: Quiescent phase. The SED 
of AG Dra during quiescent phase is characterized by (i) a dom- 
inant giant's radiation in the optical, which results in large dif- 
ferences between amplitudes of the orbitally-related variation 
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Fig. 13. The SED of AGDra during quiescence (left panels) and active phases (right). Note the variation of the nebular emission and, in 
contrast, the stability of the hot stellar component with the orbital phase during quiescence. The SED during the active phases is characterized 
by a significant increase of the nebular emission, which constrains very high temperature of the ionizing source. Dates of the IUE observations 
are marked in the U/V LCs (top). 



(AU » AB ^ AV), (ii) a strong nebula radiating at a high tem- 
perature (T e 20 000 K) and varying in its EM along the orbit, 
and (iii) a stable hot stellar component of radiation, which is 
contrary to the eclipsing systems. Fitting the spectrum from 
the maximum (27/10/89) we found the lower limit of the hot 
star temperature, T™ n = 110 000 K (k h = 5.4 10~ 25 ), to sat- 
isfy the ionization condition given here through the parameter 
(5=1 (Sect. 4. 1, Eq. 23). These parameters determine the lower 



limit of the hot star's ionizing capacity to produce the observed 
emission measure in maximum given by the fitting parameters 
k N = 5.5 10 14 ctrT 5 and T e = 21 500 K. 

Radiation from the ultraviolet: Active phase. The SED dur- 
ing the major 1994 outburst is dominated by the nebular radi- 
ation at a high electron temperature, T e = 35 000 + 5 000 K 
(£n = 2.1 + 0.4 10 16 ctrT 5 ) which determines the continuum 
profile in the ultraviolet and significantly affects the optical 
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strong and variable nebular emission dominates the ultraviolet while 
in the optical/IR we can detect radiation from only the red giant. 
Variation in the U band reflects probably that in the nebular radiation. 
Optical photometry is from Hric et al. (1994). 



with a large contribution in the infrared. At 6 cm radio wave- 
lengths, Ogley et al. (2002) measured a total emission of 
~ 1 mJy from AG Dra in 2000 March. They suggested a neb- 
ular nature for this emission. Also during the smaller eruptions 
(e.g. 1985, 1986) the nebular emission represents a significant 
component of radiation in the UV spectrum (right panels of 
Fig. 13). Here we present the case of the 1985 eruption. We 
estimated the lower limits for the temperature of the ionizing 
source, T™° = 150 000 K and 180 000 K, for the major 1994 
'cool' outburst and the 'hot' 1985 eruption, respectively. 



4.14. CQDraconis 

Reimers et al. (1988) discovered a UV bright companion of 
the M3III giant 4 Dra, determined elements of the spectro- 
scopic orbit of the giant and suggested that the variation in the 
UV could be produced by a CV of an AM Her- type. Recently, 
Wheatley et al. (2003), based on their X-ray observation, sug- 
gested that CQ Dra is most likely a symbiotic binary contain- 
ing a white dwarf accreting material from the wind of the 
red giant. Here we selected IUE observations from two dates 
(SWP24265 + LWP04630, 23/10/84, tp = 0.80 and SWP29421 
+ LWP09295, 11/10/86, y> = 0.22) to verify if the observed 
properties of the UV continuum are consistent with the ion- 
ization model of symbiotic binaries. 

Radiation from the giant. The red giant in CQ Dra is a 
dominant source of radiation in the optical/IR wavelengths 
(Fig. 14). Therefore we used also average values of the UBV 
measurements (U = 8.35, B = 6.6 and V = 5.0 mag) together 
with the JHK photometry of Kamath & Ashok (1999) to se- 
lect an appropriate synthetic spectrum. In addition, we tried 
to match an 'emission' bump from the giant at ~ /13 200A 
observed by the IUE. In this manner we selected the syn- 
thetic spectrum of T e g = 3 700K and determined the scal- 
ing factor kg = 1.5 10~ 16 , which yields the observed bolomet- 
ric flux F° hs = 1.6 10 6 erg cm" 2 s _1 . For the HIPPARCOS 
distance of 178+15 pc (Perryman et al. 1997) the giant's ra- 
dius R g = d x 9 g = 95 (of/ 178 pc)/?© and its luminosity, 
Lg = 1600 (d/ 178 pc) 2 L . 

Radiation from the ultraviolet. The profile of the ultravio- 
let continuum with the superposed emission features is typical 
of a symbiotic binary during quiescent phase. The continuum 
is nearly flat between approximately 1 600 to 3200 A, which 
suggests a dominant contribution from the nebula radiating at a 
high electron temperature. This constrains a high temperature 
for the hot stellar radiation to produce the observed nebular 
emission. In addition, the two selected spectra show a varia- 
tion in the emission measure by a factor larger than 2, while 
the scaling of the hot stellar radiation persists practically un- 
changed. This implies an increase of the hot star temperature 
from about 7/™" = 110 000 K to about 145 000 K. Parameters 
of our solution (a low luminosity, = 6.6 - 14 (J/178 pc) 2 L© 
and the effective radius, Rf = 0.007 - 0.006 (d/ 178 pc)/?©) 
suggest that the sole source of the energy detected from the 
ultraviolet is accretion from the giant's wind at the rate of 
M acc = a few x lO~ 9 M0yr _1 for the mass of a white dwarf, 
M WD = 1.0 M© (Hamada & Salpeter 1961). According to de- 
scription on wind accretion onto white dwarfs (Livio & Warner 
1984) these accretion rates require the mass-loss rate from 
the giant Mw 10~ 7 MQyr~' for the separation of the bi- 
nary components at times of observations, 750 and 900/?© 
(Mj = 2.5 M© and elements from Reimers et al. (1988)). Then 
an increase in Mw gives a larger M acc and a higher temperature 
of the ionizing source (inner parts of an accretion disk), which 
causes an increase in the emission measure, because a larger 
fraction of the giant's wind can be ionized. The first results of 
studying the accretion process in CQ Dra are discussed in more 
detail by Skopal (2005). 
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Fig. 15. The SED of V443 Her at two different positions of the binary 
(marked in the LC top). The nebular component of radiation varies 
with the orbital phase, while the hot stellar radiation is stable. 



Radiation from the giant. The JKLM photometry 
(Kolotilov & Yudin 1994) and a flat optical BVR continuum 
suggest a low effective temperature for the giant's photosphere 
(Fig. 15). A synthetic spectrum of L eff = 3 300K scaled to 
F° bs = 1.19 lO^ergctrrV 1 (k g = 1.8 1(T 18 , 6 g = 1.33 1(T 9 ) 
satisfies best the photometric measurements. Also the surface 
brightness relation for M-giants gives a very close value of 6 g = 
1 .37 1CT 9 for the reddening-free magnitudes, K = 5.32 and J = 
6.56 mag. As there is no reliable estimate of the giant's radius in 
the literature we adopted R g ~ 130/?q according to the empir- 
ical relations between the radius and T e g for M-giants (Belle et 
al. 1999). Then the distance d = R g /8 g ~ 2.2(/? g /13Ofl )kpc 
and the luminosity L g ~ 1 800 (c//2.2kpc) 2 Lq. 

Radiation from the ultraviolet. The stellar component of 
the radiation represents a dominant contribution in this spec- 
tral region. It has a relatively small slope in the far-UV and 
rivals other components of radiation up to the B band (Fig. 15). 
Our model of the 21/10/80 spectrum suggests a low temper- 
ature, T h = 39 000 + 11000/ - 5 000K, where uncertain- 
ties were determined in the same way as in Skopal (2001b). 
Both the spectra require the same scaling, = 2.32 10 -23 
(#h = 4.82 10~ 12 ), which corresponds to the effective ra- 
dius Rf = 0.47(if/2.2kpc)fl© and the luminosity L h = 
460(ii/2.2kpc) 2 Lq. On the other hand, the nebular emission 
is faint and varies along the orbit: Ic^ = 1.0 and 0.4 10 15 ctrT 5 
at tp = 0.5 (r e = 15 500 + 3 000K) and <p = 0.25 (T e = 
12 000 + 2 000K), respectively. However, the corresponding 
parameter S = 3.9 and 2.0 is not consistent with the fitting pa- 
rameters Th, kh and Ic^. The observed hot stellar radiation is 
not capable of producing the nebular emission, because of the 
too low temperature given by the fit (Lh «: r™ n ~ 75 000). 
This conflicting situation signals that an unseen part of the ion- 
izing source has to radiate at a higher temperature to produce 
a surplus of the nebular flux in addition to that generated only 
by its observed portion. Qualitatively, this could be understand- 
able if the accreted matter has a disk-like structure, whose outer 
rim occults the central hot core in the direction to the observer. 
However, this would suggest a rather high inclination of the 
V443Her binary orbit. Parameters Lh and R^ derived from 
our models represent the lower and upper limits of their real 
values. Therefore in Table 3 we also present a formal solution 
for L h = r™ n = 75 000 K (Sect. 4.1, Eq. 23). 



4.15. V443Herculis 

V443 Her is a quiet symbiotic star - no outburst has been 
recorded to date. The LC displays a marked wave-like varia- 
tion of the optical light along the orbital phase (Fig. 15). Recent 
studies of this symbiotic (Dobrzycka et al. 1993; Kolotilov 
et al. 1995; Skopal 1996) suggested a low orbital inclina- 
tion (« 30°, and 18°, respectively). We selected two pairs of 
IUE spectra taken at the maximum of the light (SWP 10439 
+ LWR09118, 20/10/80, <p = 0.50) and at the quadrature 
(SWP42772 + LWP21540, 20/10/91, y = 0.25). No well- 
exposed spectrum at/around the minimum was available. 



4.16. YYHerculis 

Munari et al. (1997b) demonstrated that the LC of YY Her 
underwent four main outbursts and a few bright stages. 
Photometric measurements made after the main 1993 outburst 
(Hric et al. 2001; Mikolajewska et al. 2002a) revealed a 
relatively deeper and narrower primary minima than prior to 
the outburst (Fig. 16), which suggests a high inclination of 
the orbital plane. In addition, a secondary minimum in the 
VRI LCs appeared (Hric et al. 2001). Such evolution is simi- 
lar to that of AX Per during the transition from its 1989 out- 
burst to quiescence (Fig. 3 of Skopal et al. 2001a, Fig. 21 
here). We selected IUE spectra from the maximum of the 1993 
outburst (SWP47861 + LWP25733(4), 13/06/93, (p = 0.37), 
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Fig. 16. The SED of YY Her during quiescence and two different 
phases of activity. Visual LC is composed from the data of Munari 
et al. (1997a) (o) and those available from the CDS database (•). The 
eclipse effect observed in 1994.43 is shown in detail. Crosses (x) in 
the SED represent continuum flux-points from the optical spectrum 
taken on 10/06/93 (Munari et al. 1997a). 



a bright phase in 1981 (SWP15652 + LWR12080, 05/12/81, 
<p = 0.27) and quiescence (SWP09773 + LWR08493, 12/08/80, 
<p = 0.46). 



Radiation from the giant. Relatively faint RIJKL magni- 
tudes of a normal giant (e.g. Miirset & Schmid 1999) suggest 
that YYHer is a distant object. The corresponding synthetic 
spectrum is characterized by T eS = 3 500 K and 8 g = 3.9 10~ 10 
(i.e. F° bs = 1.28 10~ 9 erg cm -2 s -1 ). The average reddening- 
free values of 7=8.98 and K-7.S9 provide nearly identical 
angular radius, 6 g = 4.0 10~ 10 (Dumm & Schild 1998). We es- 
timated the radius of the giant from the profile of the 1994 min- 
imum in the smoothed visual LC (Fig. 16). This can be inter- 
preted in terms of the eclipse of the hot star pseudophotosphere 
by the giant. Its timing (f\ ~ JD 2 449486, t 2 ~ JD 2 449492, 
h ~ JD 2 449522, f 4 ~ JD 2 449537 : ±2-3 days) determines 
the linear size of the giant's stellar disk that eclipses the ob- 
ject, /?g = 0.21 A. This parameter can be expressed through the 
radius of the giant, R g , and the orbital inclination ;' as 

^ g /A) 2 -cos 2 (0. (25) 



R E JA 



With the aid of the empirical relations between the radius 
and effective temperature for giant stars (Belle et al. 1999) 
we adopted R g ~ 110 Rq, which satisfies Eq. (25) for i = 
80° and A = 403 Rq (M t = 2.5 Mq). This yields the dis- 
tance d = 6.3 (/?g/110/?©)kpc and the giant's luminosity L g = 
1 600 (d/6.3 kpc) 2 Lq. On the other hand, based on the pres- 
ence of a secondary minimum in the VRI LCs, Mikolajewska et 
al. (2002a) suggested a Roche-lobe filling giant in the system, 
which corresponds to the distance d = 10 ± 3 kpc. However, 
the / LC shows rather narrow minima (see Fig. 2 of Hric et al. 
2001) than can be ascribed to the wave-like variation due to 
the tidal distortion of the giant. Also a statistical approach of 
Miirset & Schmid (1999) suggests the giant in YYHer to be 
well inside its Roche lobe. 

Radiation from the ultraviolet: Quiescent phase. Radiation 
of YY Her during quiescence comprises a hot stellar compo- 
nent of radiation of r n > 105 000 K = T™ 11 , the angular radius 
0h < 4.9 10~ 13 and the nebular component of radiation char- 
acterized by T e = 21 500 + 3 000K and k N = 2.4 10 14 ctrT 5 . 
Other derived parameters are given in Table 3. 

Radiation from the ultraviolet: Active phases. YY Her is the 
only case in our sample of objects that shows a very different 
continuum profile during active phases. During the maximum 
of the star's brightness (1993) the profile of the SED is of the 
same type as observed for AGDra (cf. Fig. 13), while during 
the 1981 bright phase it was similar to that of all other active 
symbiotics with a high orbital inclination (e.g. CICyg). In the 
former case the SED from the ultraviolet to the near-IR region 
is dominated by the nebular radiation at a high electron tem- 
perature (T e ~ 40 000 K, k N = 2.1 10 15 ctrT 5 ), which is pow- 
ered by a strong stellar source radiating at Th > 160 000 K = 
T™ n and scaled with <9 h < 6.3 10~ 13 (i.e. 6 < 1). The profile 
of the simultaneously observed optical spectrum (Munari et al. 
1997a) agrees well with the modeled SED (Fig. 16). In the lat- 
ter case the UV continuum results from superposition of radi- 
ation from a cool stellar pseudophotosphere (Th = 22 000 K, 
0h = 3.2 10~ 12 ) attenuated with n H ~ 1-4 10 23 crrT 2 , a low- 
temperature nebula (r e = 13 000 K, fe N = 1.5 10 14 crrT 5 ) and a 
high-temperature nebula (T e = 40 000 K, k N = 1.3 10 14 ctrT 5 ). 
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Fig. 17. Top: Visual and U LCs of RW Hya (Skopal et al. 2002b) with 
added (/-magnitudes estimated from available IUE observations. The 
arrow marks position of the used spectrum. Middle and bottom: The 
SED near to the optical maximum shows typical characteristics of a 
quiescent phase with a high luminosity of the hot object, while during 
the eclipse only a very faint high-temperature nebular emission can be 
recognized. 



4.17. RWHydrae 

RW Hya is a stable symbiotic system for which no optical erup- 
tion has been observed. Its LC displays a wave-like variation 
as a function of the orbital phase (Fig. 17). A high inclina- 
tion of the orbital plane is indicated by the eclipse effect due 
to attenuation of the far-UV continuum by Rayleigh scattering 
(e.g. Schild et al. 1996). IUE spectra of RWHya have been 
studied in detail by many authors (Kenyon & Mikolajewska 
1995; Schild etal. 1996; Dumm et al. 1999; Sion et al. 2002). 



We selected one observation around the maximum (SWP05863 
+ LWR05110, 19/07/79, cp = 0.32) and one around the mini- 
mum of the optical light (SWP31039 + LWP10825, 24/05/87, 
ip = 0.06) to demonstrate our approach for this case. 

Radiation from the giant. A synthetic spectrum of r e jf = 
3 800K with scaling of k g = 2.6 10~ 18 (<9 g = 1.61 10 9 , 
F° bs = 3.10 10 s erg cm" 2 s" 1 ) matches well the flux-points 
corresponding to the BVRIJHKL photometry. Reddening-free 
magnitudes, K = 4.65 and J = 5.72 mag, also give a similar 
value of 6» g = 1.77 10~ 9 (Dumm & Schild 1998). The giant 
dominates the spectrum from the B band to the infrared. This 
results in a large amplitude difference of the wave-like vari- 
ation at different passbands (AU » AV, Fig. 17). Thus, the 
nebular emission, which is responsible for this type of variabil- 
ity, is not able to rival the giant's contribution from B to longer 
wavelengths. Assuming the synchronous rotation of the giant 
with the orbital revolution, Schild et al. (1996) determined its 
radius to R g = 58.5 + &Rq. This quantity then yields the dis- 
tance d = R g /6 g = 820 + 112(fl g /58.5/?0)pc and the lumi- 
nosity, L g = 650 + 180 (cZ/0.82 kpc) 2 Lq, where uncertainty in 
only R g is included. These parameters are similar to those de- 
termined by Schild et al. (1996). 

Radiation from the ultraviolet. A steep far-UV continuum 
(practically within the SWP spectrum) reflects a dominant con- 
tribution from the hot stellar source in this region. On the other 
hand, a flat near-UV continuum of the LWR spectrum signals 
a significant nebular radiation here - an illustrative example 
of a quiescent symbiotic's UV spectrum. Our best fit of the 
19/07/79 spectrum corresponds to T h = 110000 + 30 000 K, 
k h = 6.3 10~ 24 and T e = 15 000 + 2 000K with scaling, 
&n = 2.8 10 15 ctrT 5 . This solution gives the parameter 6 = 0.54, 
which means that a fraction of ionizing photons is not con- 
verted to the nebular radiation. In the sense of the STB model 
the ionized zone is open. The situation in which all the Lph 
photons are involved in the ionization process corresponds to 
T min _ go 000 K, which represents the lower limit of the hot 
source temperature in RW Hya. It is also close to the Zanstra 
temperature of 75 000 K derived by MNSV The spectra of 
RW Hya taken close to the giant's inferior conjunction suffer an 
additional (to the Rayleigh scattering) wavelength-independent 
attenuation (Dumm et al. 1999) as observed for other symbi- 
otics with a high orbital inclination. During the totality only a 
high-temperature nebula can be recognized. 

4.18. SYMuscae 

SY Mus is another stable system with a high luminosity of the 
hot component. Its visual LC displays a strictly periodic wave- 
like variation with the orbital phase (e.g. Pereira et al. 1995). 
Figure 18 also shows a variation in the emission measure along 
the orbital phase scaled to the effective wavelength of the U 
band as suggested by Skopal (2001a). The emission measure 
was estimated from the near-UV, where the nebula represents 
a dominant source of radiation (Eq. 18, T e = 18 500 K). The 
system is eclipsing. When the binary approaches the position 
of the inferior giant's conjunction, the Rayleigh scattering sig- 
nificantly affects the far-UV radiation (e.g. Pereira et al. 1995). 
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Fig. 18. Top: Visual and U LC of SYMus. The former was compiled 
from visual estimates available at the CDS database. The latter was 
reconstructed from the IUE low-resolution spectra by the same way 
as suggested by Skopal (2001a). Middle and bottom: Observed and 
reconstructed SEDs at the position near to its optical maximum. 



Near to the conjunction (<p — + 



practically flat and at the level of 1 - 2 10 13 erg cm 2 s 



0.05) the UV continuum is 

'A-'. 

Also in this case an additional wavelength-independent con- 
tinuum attenuation can be recognized (Table 2 of Dumm et 
al. 1999). As a result we selected only one observation ex- 
posed at the position near to the maximum, where the effect 
of the additional opacity source can be assumed to be minimal 
(SWP15594 + LWR12059, 29/11/81, y = 0.62). 

Radiation from the giant. In good agreement with the spec- 
tral classification of the giant in SYMus (M4.5III, Miirset 
& Schmid 1999) a synthetic spectrum of r eff = 3 400K 
matches well photometric measurements in the IR. Its scaling, 
k g = 3.6 10~ 18 , (0 g = 1.9 10~ 9 ) corresponds to the bolomet- 
ric flux, F° bs = 2.6 10~ 8 erg ctrT 2 s _1 . According to the surface 
brightness relation for M-giants (Dumm & Schild 1998), the 
reddening-free magnitudes, K = 4.55 and J = 5.75 mag, give 
the same value of 9 g . These fitting parameters define the dis- 
tance d = 1.0 + 0.15 (/f g /86/?o) kpc and the luminosity L g = 



850 + 250(^/1 kpc) 2 L , where the radius R g = 86 + 13 Rq 
was determined by Schmutz et al. (1994) in the same way as 
forRWHya. 

Radiation from the ultraviolet. There is a striking similarity 
in the SEDs of SYMus and RWHya (Figs. 17 and 18), which 
suggests very similar ionization conditions in these systems. 
The evolution in the Ha line profile along the orbital motion 
is qualitatively of the same type in both systems (Schmutz et 
al. 1994; Schild et al. 1996). Our solution for the selected 
spectrum corresponds to T h = 110000K (= Zanstra temper- 
ature), k h = 4.8 10~ 24 and T e = 18 500 + 2500K with scal- 
ing, fc N = 2.5 10 15 cm" 5 , which gives the parameter 6 = 0.55. 
These parameters can be discussed in the same way as those 
for RW Hya. Here the lower limit of the hot source temperature 
is T™ n = 80 000 K. The only larger difference is in the rela- 
tive contributions from the giant and those produced by the hot 
star in the U band, which is due to the difference in T e g. This 
situation makes the amplitude A£/(SYMus) > A£/(RWHya) 
despite that orbital inclinations of both systems are compara- 
ble. As a by-product of fitting the UV spectra, we determined 
a new value of reddening to SYMus, E B -v = 0.35. A larger 
quantity of £b-v = 0.5 - 0.45, previously suggested by MNS V 
and Pereira et al. (1995), yields too steep a far-UV continuum 
(from « 1 350 A), which is not possible to match by any rea- 
sonable black-body radiation. The 2 100 - 2 300 A spectral re- 
gion displays a higher continuum level than that at surrounding 
wavelengths if one consider the influence of the iron curtain. 
We note that this region is rich in strong absorption features 
(see in detail the well-exposed spectra, e.g. LWR12059). 



4.19. ARPavonis 

Skopal et al. (2001b) recently summarized the historical, 1889 
- 2001, LC of ARPav. It is characterized by about 2 mag deep 
minima - eclipses - and strong out-of-eclipse variations of 
between about 12 and 10 mag. This behaviour suggests that 
AR Pav persists in an active phase. The UV continuum shows 
the same type of profile during the whole IUE mission (e.g. 
Fig 1 of Skopal 2003b). Therefore we selected observations 
well outside the giant's inferior conjunction (SWP13956 + 
LWR10570, 10/05/1981, <p = 0.37) and just in the eclipse 
(SWP16949 + LWR13235, 13/05/1982, <p = 0.98). 

Radiation from the giant. We compared the RIJHKL fluxes 
with a synthetic spectrum of r e tf = 3 400K and 9 g = (6.4 ± 
0.5) 10~ 10 , which corresponds to the bolometric flux F° bs = 
(3.1 + 0.3) 10~ 9 erg ctrT 2 s _1 . The average reddening-free val- 
ues of J = 7.96 and K = 6.86 mag yield 9 g = 6.5 10~ 10 (Dumm 
& Schild 1998). Schild et al. (2001) derived the distance d = 
4.9 kpc from the 'rotational' giant's radius, R g = (130 ± 25) Rq 
and the brightness surface relation for M-giants. This value is 
in particular agreement with that derived from a 106-day pe- 
riod of pulsations of the giant seen in some parts of the vi- 
sual LC (Skopal et al. 2000). Our value of 8 g then gives the 
radius of the giant R g = (139 + 10) Rq, which agrees with 
that determined by Schild et al. (2001). The uncertainty results 
only from that in 6 g . This quantity corresponds to the giant's 
luminosity, L g = (2 300 + 400)(t//4.9kpc) 2 L Q - According to 
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Fig. 19. The SED for AR Pav out and in the eclipse (positions of the 
spectra used are marked by arrows in the relevant part of the LC). 
Its out-of-eclipse profile in the ultraviolet resembles that observed in 
other eclipsing systems during active phases. 



the orbital solution, the separation between the components, 
Asin(z') = 457 + 35% (Quiroga et al. 2002) and the lin- 
ear size of the giant's stellar disk that eclipses the object, 
Rg/A = 0.30 + 0.02 (Skopal et al. 2000) limit the R g radius 
to between 207 + 25/?© (R g ~ R h , i = 72°) and 137±20fl© 
(i = 90°) (Eq. 25, see also Quiroga et al. 2002). From this 
point of view, R g = (139 ± 10) Rq corresponds to the orbital 
inclination i = 87 + 3°. 

Radiation from the ultraviolet. The profile of the UV 
continuum on the out-of-eclipse spectra is of the same type 
as observed for other active symbiotics with a high orbital 
inclination: A Rayleigh attenuated far-UV continuum (here 



n H ~ 4.5 10 22 cm 2 ), a rather cool radiation from the hot 
stellar source (T b ~ 22 000 K, k h = 7.9 10~ 23 ), which is 
markedly affected by the iron curtain absorptions, and the 
nebular components of radiation from both the LTN (r e = 
15 500 K) and HTN (T e ~ 40 000 K). Radiation from the 
HTN is identified in all spectra. During eclipses, as the only 
light visible in the ultraviolet and on the out-of-eclipse spec- 
tra it is indicated by a non-zero level of the Rayleigh atten- 
uated continuum near to the Lya line. Its level is about 1 - 
2 x 10~ 13 erg crrT 2 s _1 A -1 , which corresponds to a large emis- 
sion measure, EMutn = 2.9 10 60 (<f/4.9kpc) 2 cirr 3 and/or the 
luminosity Lhtn = 530 (d/A.9 kpc) 2 Lq. Our solution is plotted 
in Fig. 19. 



4.20. AGPegasi 

AG Peg is the slowest symbiotic nova. In the mid- 1850 it rose 
in brightness from ~ 9 mag to ~ 6 mag and afterwards followed 
a gradual decline to the present brightness of ~ 8.7 mag in V. 
Evolution in its spectrum during the last century was described 
by many authors (e.g. Kenyon et al. 1993). The orbitally re- 
lated wave-like variation in the LC developed around 1940 
(Fig. 6 of Skopal 1998, and references therein). From that time 
to the present AG Peg faded from about 7.5 to 10 mag in the 
B band. Therefore we selected two earlier observations within 
one orbital cycle, but at/around different conjunctions of the 
components (SWP03830 + LWR03376, 05/01/79, <p = 0.96; 
SWP07407 + LWR06390, 14/12/79, (p = 0.38) to show the 
orbitally-dependent nebular emission at a constant hot stel- 
lar radiation, and one observation taken later when the star's 
brightness significantly faded (Y1JO0308T + Y1JO0403T + 
Y1JO0309T + Y1JO0406T, 13/11/93, (p = 0.63). This was 
made by the Faint Object Spectrograph (FOS) on the board the 
Hubble Space Telescope. 

Radiation from the giant. The cool component in AG Peg 
was classified as a normal M3 III red giant (Miirset & Schmid 
1999; Kenyon et al. 1993). Its rather bright IR magnitudes sug- 
gest it cannot be very distant. The flux points corresponding 
to the RIJKLM measurements can be matched by a synthetic 
spectrum with r eff = 3 600 K, log(g) = 1.0 and 6 g = 2.4 10~ 9 , 
which yields the flux F° bs = 5.6 10~ 8 ergcm^s -1 . Our mean 
dereddened magnitudes, K = 3.85 and J = 4.98 mag, also give 
a similar value of 9 g = 2.6 10~ 9 . Kenyon et al. (1993) deter- 
mined the distance modulus for bolometric magnitudes under 
the assumption that the giant in AG Peg is similar to other red 
giants. They obtained d = 800 pc, which corresponds to R g = 
dx0 g = 85RQ and its luminosity L g = 1 100 (c//800pc) 2 Lq. 
These quantities are identical to those determined by Kenyon et 
al. (1993). Other estimates in the literature also rely on com- 
mon properties of red giants as a class. For example, MNSV 
employed the same method, but for K magnitudes, and ob- 
tained d = 650 pc. As there is no independent determination 
of the giant's radius in AG Peg we adopted that published by 
Kenyon et al. (1993), which conforms better to our solution. 

Radiation from the ultraviolet. Properties of the UV spec- 
trum and optical LCs during the investigated period resem- 
ble those of classical symbiotic stars during quiescent phases 
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(Fig. 20). In addition, the overall fading of the star's brightness 
results from that in (see Fig. 4 of Miirset & Nussbaumer 
1994, and Fig. 20 here). According to Miirset & Nussbaumer 
(1994) a constant hot star temperature T h = 95 000 K can be as- 
sumed during the period 1978 - 1993. We adopted this quantity 
in our modeling the SED. The following points are relevant. 

(i) The 1979 spectra show the phase-dependence of the neb- 
ular radiation, which follows that in the LCs. The observed EM 
increased by a factor of about 2 from 05/01/79 (tp = 0.96) to 
14/12/79 (<p = 0.38). In contrast to, e.g. RW Hya and SYMus, 
the far-UV continuum is not subject to orbitally-related varia- 
tion, which suggests that the AG Peg orbit is not highly inclined 
to the observer. 

(ii) We estimated the lower limit of the hot star temperature, 
r min = 70 000 K for the 14/12/79 spectrum. The H II region in 
AG Peg is very open (Fig. 11 of Kenyon et al. 1993), which 
means that a significant fraction of the L p h photons escapes the 
system. This implies Th » T™ n to produce the observed EM 
for an appropriate fit. For example, if only one half of the ion- 
izing photons gives rise to the observed nebular emission then 
the hot star temperature has to increase to Th = 100 000 K to fit 
the far-UV continuum (Eq. 23 for S = 0.5). 

(iii) The situation for AG Peg is complicated by its hot star 
wind (Voo -900 - 1000 km s" 1 , M ~ 3 10 _7 M G yr-\ Vogel 
& Nussbaumer (1994)), which represents an additional source 
of emitters contributing to the observed nebula. However, the 
emission measure from the hot star wind, £Mw, is negligi- 
ble with respect to that observed. We estimated its quantity to 
be EM W = 8.3 1 58 cm" 3 by using relation (8) of Skopal et al. 
(2002a) for/?* = Rf = 0.18%, the radius r min ~ 1.5 R*, 
from which the wind becomes optically thin (Leitherer 1988) 
and the parameter /J = 1 as assumed for hot star winds. Thus 
£M W <K EM° hs . This suggests that a major part of the nebular 
emission has to come from a more dense region in the system, 
which could be associated with the colliding region of the two 
winds located probably at the vicinity of the giant's hemisphere 
facing the hot star (Kenyon et al. 1993). 

(iv) Finally, the 1993-spectrum shows a higher efficiency 
in producing the nebular component of radiation (6 = 0.73). A 
decrease of the mass-loss rate from the hot star by a factor of 
~ 5 (Vogel & Nussbaumer 1994) results in a larger expansion 
of the massive giant's wind towards the hot star and thus in 
a larger opening of the colliding zone. This makes it possible 
for a larger fraction of the L p h photons to be involved in the 
ionization/recombination process. A quantitative model should 
support this view. 

4.21. AXPersei 

The eclipsing nature of this classical symbiotic star was unam- 
biguously revealed by photometric observations of the ~ 1 mag 
deep narrow minimum - eclipse - in the UBV LCs during 
its recent 1989 major outburst (Skopal 1991). The historical 
LC of AX Per is characterized by long-lasting periods of qui- 
escence with superposition of a few bright stages (Skopal et 
al. 2001a). Accordingly we selected four IUE spectra taken 
at very different brightness/activity levels: During quiescence 
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Fig. 20. The SED for AG Peg for three selected dates at different 
brightness of the system. Positions of the used spectra are marked by 
arrows in the LC (top). Data are from Luthardt (1984, 1989), Komarek 
(1990), Belyakina (1992), Tomov & Tomova (1998, 2001), Skopal et 
al. (2004) and from 2003 our unpublished observations. 



(SWP24278 + LWP04619, 23/10/84, ip = 0.71), during the 
1978 outburst (SWP03755 + LWR03332, 31/12/78, <p = 0.59), 
during the 1990 eclipse (SWP40238 + LWP19327, 30/11/90, 
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Fig. 21. The SED for AX Per during the activity, quiescence and 
eclipse. Positions of the spectra used are marked by arrows in the LC. 
Small dots represent means of the CDS visual data in the 30-day bins. 



ip = 0.99) and during the transition to quiescence (SWP49127 
+ LWP26699, 07/1 1/93, y = 0.57). 

Radiation from the giant. Infrared photometry of AX Per 
was summarized by Skopal (2000), who derived the follow- 
ing parameters of the giant in the same way as in this pa- 
per: T eB = 3400 ± 150 K, 6» g = (1.3 + 0.2) 10" 9 , F° hs = 
(1.4 + 0.3) 10~ 8 erg cm -2 s -1 . The angular radius then gives the 
distance d = 1 700 + 200 (RJ 102 Rq) pc, where the giant's ra- 
dius was derived from the eclipse timing under the assumption 
that ;' = 90°, which is supported by a rectangular profile of the 
minimum (Skopal 1994). The giant's corresponding luminos- 
ity L g = 1200 + 4OO(c//1.7kpc) 2 L . 

Radiation from the ultraviolet: Quiescent phase. The best 
solution of the 23/10/84 spectrum yields a low temperature of 
the observed hot stellar source, T h = 66 000 + 15 000/ -7 000 K 
(k h = 1.7 10~ 24 ), but a strong nebular emission characterized 
by T e = 25 000 ± 3 000 K and scaled with k N = 7.2 10 14 ctrT 5 . 
These parameters give 6 = 2.3 (Eq. 21), i.e. the corresponding 
amount of the L p h photons is not capable of producing the ob- 
served quantity of the nebular emission. This situation signals 
the presence of a material around the accretor of a disk-like 
structure, which makes it optically thick to a larger distance 
from the hot surface of the central star (> /?jj ff ) in directions on 
the orbital plane (i.e. on the lines of sight; i ~ 90°). This means 
that in AX Per the accreted matter has a disk-like structure even 
during quiescent phases. For comparison, we also present a for- 
mal solution for 6=1, i.e. T h = 7J in =100000K (Table 3). 
The solution for this case gives = 0.065(d / U kpc) Rq 
and L h = 380(J/1.7 kpc) 2 LQ, which are similar to those deter- 
mined by MNS V. Our somewhat larger value of Lh is probably 
caused by different scaling, because we considered the influ- 
ence of the iron curtain absorptions. 

Radiation from the ultraviolet: Active phases. It is charac- 
terized by a flat profile with signatures of the Rayleigh atten- 
uated far-UV continuum and the influence of the iron curtain 
absorptions. An appropriate discussion of such a spectrum was 
introduced, for example, in Sect. 4. 10 for CI Cyg. Also here the 
three-component model, composed of the HSS, LTN and HTN 
radiative contributions, has to be applied to get a satisfactory 
fit to observations. Fitting and derived parameters of our solu- 
tions for both selected spectra (31/12/78, 07/11/93) are given 
in Table 4 and Fig. 21. 

Radiation from the ultraviolet: Eclipse. The UV spectrum 
during the 1990-eclipse can be matched by the HTN radiating 
at T e = 32500 + 3 000K and scaled with k N = 8.0 10 14 ctrT 5 . 
The same component was used in fitting the 31/12/78 spectrum 
to fill in the Rayleigh attenuated far-UV continuum. However, 
during the 1994-eclipse (SWP52027 + LWP29094, 04/09/94, 
(p = 0.01), at transition to quiescence, the level of the HTN 
emission was a factor of about 2.5 fainter. This suggests that the 
amount of the HTN contribution is a function of the activity. 



4.22. BD-21°3873 

Smith et al. (1997) presented the largest set of uvby photometry 
of this symbiotic star. The LCs show a pronounced orbitally- 
related variation with superposed double-sine patterns in y, b 
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Fig. 22. Approximate SED for BD-21°3873 given by its SWP spec- 
trum in the ultraviolet and UBVRIJK fluxes from photometric mea- 
surements (•). The value of U(tp = 0.31) from Munari et al. (1992) 
was extrapolated to U(tp = 0.5) according to the profile of the u- 
LC. Crosses (x) represents flux-points from the spectrum taken on 
24/02/88 (<p = 0.61) by Niehues et al. (1992). 



and v. There are no signatures of an active phase (also Fig. 9 of 
Skopal et al. 2002b). As for LTDel only SWP spectrum was 
available in the IUE archive (SWP40919, 26/02/91, (p = 0.49). 
We determined a new value of Lb-v = 0.20 by the same way 
as for LTDel (Sect. 4.12). 

Radiation from the giant. Fluxes of the BVRIJK photomet- 
ric measurements can be compared well with a synthetic spec- 
trum calculated for T e g = 4 300K with agreement with the 
atmospheric analysis of Pereira & Mello (1997) and Smith et 
al. (1997). The scaling, k g = 2.4 10~ 19 , then gives the bolomet- 
ric flux F° bs = 4.7 10~ 9 erg cirT 2 s _1 . From the surface gravity, 
log(g) = 1.69 ± 0.17 (Pereira & Mello 1997), one can derive 
the giant's radius as R g = (29 ± 6)(M g /1.5Mo) 1/2 fl and the 
distance d = R s /8 g = (1.3 ± 0.3)(M g /1.5M Q ) 1/2 kpc, where 
the giant's mass was adopted as for AGDra, because of their 
very similar atmospheric parameters (Smith et al. 1997). The 
giant's luminosity L g = 260 (d/1.3 kpc) 2 Lq. 

Radiation from the ultraviolet. Observed and approximate 
modeled SEDs are shown in Fig. 22. As for LT Del we adopted 
T h = 1 00 000 K, which requires a scaling of k h = 1.71 0~ 24 and 
suggests parameters for the nebular emission, T e ~ 15 000K 
and &n ~ 3.3 10 14 cirT 5 . This solution gives 6 = 0.36. Also 
T h = 75 000 K (Schmid & Nussbaumer 1993) is above T™ n , 
because the corresponding 6 = 0.63. This results suggest an 
open ionized zone (X > 1), which is important in calculating 
the mass-loss rate from the giant (Sect. 5.2.1). 

5. Summary 

By disentangling the composite spectrum of 21 S-type sym- 
biotic stars we isolated four basic components of radiation 
contributing to their UV/optical/lR continuum: Radiation from 
the (i) cool giant, (ii) hot stellar source (HSS), (iii) low- 
temperature nebula (LTN) and (iv) high-temperature nebula 
(HTN). Physical parameters of these components are intro- 
duced in Tables 2-4 and corresponding solutions are shown 



Table 2. Fitting and derived parameters: Cool components 
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Useful conversions according to the section 3.1: k g = G 2 g = F^/crT^, 
R g =dx6 g , L g = 4nd 2 Ff s and d in Table 1 . 

in the relevant figures. During quiescent phases, components 
from the hot object show orbitally-related variations and unpre- 
dictable changes in T e , Lh and EM. Transition to active phases 
is followed by a dramatic change in Th and by the appearance 
of the HTN component in the spectrum. Below we summarize 
and discuss their common properties and observed effects. 

5.1. Radiation from the giant 

This component of radiation in the spectra of S-type symbi- 
otic stars is determined by two fitting parameters - the angu- 
lar radius, 6 g and the effective temperature, T ef[ (Sect. 3.1). 
Both are connected throughout the observed bolometric flux, 
^■obs _ qI^jA^ yyjjjcjj defines the giant's luminosity as 

L g = And 2 Ff\ (26) 

To test common properties of giants we examined the quan- 
tity of ,F° bs as a function of the distance d in the sense of the 
last relation. Fitting the relevant quantities from Table 1 and 2 
determines this function as 

log(F° bs ) = -2 log(rf) - (7.30 + 0.05) for red giants (27) 
and 

logCF° bs ) = -2 log(J) - (8.03 + 0.05) for yellow giants. (28) 

The flux is in erg cirT 2 s~' and the distance in kpc. The constant 
corresponds to a characteristic luminosity L g = 1 600 + 200 and 
290 + 30 Lq for red and yellow giants, respectively. Figure 23 
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plots the observed data and fits. In the figure we distinguish gi- 
ants having independently determined radii (from eclipses, co- 
rotation, atmospheric analysis, lobe-filling giant; 11 objects), 
which then provide the distance from their 6 g . For CQ Dra and 
CH Cyg we adopted their HIPPARCOS distances. This group 
of giants confirms that cool components in S-type systems are 
normal giants. Otherwise the distances were estimated from 
S T/R g dependencies (Belle et al. 1999) and the measured an- 
gular radii. CH Cyg, which is suspected to contain two red gi- 
ants in accordance with its triple-star model, was not included 
in the analysis. Finally, having the SED in the infrared, these 
relations can be used to estimate distances of S-type symbiotic 
stars. 

5.2. The quiescent phase 

5.2.1. Nebular continuum radiation 

General properties: Our solutions for 15 symbiotic stars dur- 
ing the quiescent phase showed that the nebular component 
of radiation can be characterized by a unique electron tem- 



perature. For individual objects T e runs between about 12000 
and 25 000 K with an average value of 19 000 K (Table 3). For 
Z And and AG Peg we observed a large decrease in T e by about 
10 4 K. In the case of Z And, T e decreased at approximately the 
same L N and probably also Lh, while for AG Peg this change 
was followed by a significant decrease in L N , obviously caused 
by that in Lh and thus L p h (Table 3). A decrease in T e could be 
caused by a decrease of the hot star wind, which causes removal 
of a fraction of hot electrons from the plasma created by colli- 
sions. For EG And and CQ Dra a transient enhancement of the 
emission measure observed at the same orbital phase can result 
from a relevant increase in the mass-loss rate from the giant, 
which thus supplies more emitters into the ionized zone. This 
change has a counterpart in the increase of the star's brightness 
in U. 



Mass-loss rate from the giant: The nebula during quies- 
cence arises from ionization of the cool giant wind. Skopal 
(2001a) showed that the observed emission measure is con- 
sistent with that produced by the ionization model. Therefore 



28 A. Skopal: Disentangling the composite continuum of symbiotic binaries 

Table 4. Fitting and derived parameters: Hot components during active phases 
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smaller. It represents only a fraction of the whole sphere, which 
requires a higher M w to give rise the observed quantity of EM. 
For example, mass-loss rates calculated according to Eq. (29) 
for the open H n regions characterized with X — 20 are factor 
of about 2.5 larger than those estimated from Eq. (30). Most 
of the objects we investigated here contain extended nebulae 
which is mainly due to the high number of the ionizing pho- 
tons they produce and a large separation between the binary 
components. Therefore for these objects we calculated M w ac- 
cording to Eq. (29) for X = 20. In the cases of a less opened 
ionized zone (EG And, CQDra and possibly AX Per) we in- 
tegrated Eq. (29) throughout its volume given by the parame- 
ter X. For EG And and CQ Dra, Skopal (2005) obtained M w 
and the corresponding parameter X from accretion processes 
in these binaries. For AX Per Skopal et al. (2001a) determined 
X = 1.0 from timing of the broad minimum observed during 
the 1958-72 quiescencent phase. However, as this value could 
change significantly (following the change of the LC profile), 
we present also a solution for X = 20. Results of our analy- 
sis are given in the last column of Table 3 and Fig. 24 shows 
a comparison of our values with those determined by MNSV 
and Mikolajewska et al. (2002b) from radio observations as a 
function of the hot star luminosity. The graph confirms well the 
relationship between these two parameters. Mikolajewska et al. 
(2002b) discussed this correlation in terms of the illumination 



we can calculate EM directly by integrating emission contribu- 
tions throughout the volume of the fully ionized zone as 

EM = f n 2 (r)dV, (29) 
Jh ii 

which can be determined through parameters of the giant's 
wind. The particle density n(r) satisfies the equation of continu- 
ity, M w = 4-7Tr 2 fimun(r)vw(r), in which the velocity structure 
of the wind, v w (r) = Voo(l - Rg/rf (Castor et al. 1975), v ra is 
the terminal velocity of the wind and the parameter j3 charac- 
terizes its acceleration; j3 = 2.5 satisfies a wind from a late-type 
giant (Schroder 1985). Under these conditions the upper limit 
of the EM calculated in accordance with the relation (29) can 
be expressed analytically as 

-l /m w \ 2 i r / ^r 4 i 

£M= t" j — j — 1 j 1 — j (30) 

16^m H ) 2 \ Voo / R g V \ Q) \ 

(Nussbaumer & Vogel 1989; Skopal 2001a), which is applica- 
ble for extended ionized zones in symbiotic binaries (i.e. very 
open in the sense of the STB model; parameter X » 1). This 
solution corresponds to an upper limit of the EM assuming the 
sphere around the cool star to be fully ionized from r — Q to 
r - co. The parameter Q is the location of the H i/H n boundary 
on the binary axis and the wind starts at the giant's surface. We 
comment on our calculations of Mw as follows: 

(i) The mass-loss rate estimated from relation (30) repre- 
sents its lower limit. In the real case the ionized volume is 
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heating of the outer red giant's atmoshere, which enhances the 
mass-loss rate and consequently can induce a larger Lh- 

(ii) Production of the nebular emission depends on the ter- 
minal velocity of the wind. A low-velocity wind corresponds 
to a higher concentration of emitters in the H n zone, which 
thus gives rise to the observed EM for a lower mass-loss rate. 
We adopted v M = 20 km s" 1 . Note that MNSV adopted // • Voo = 
lOkms -1 (i.e. v ro = 7kms _1 ), which provides a factor of about 
3 lower values of Mw 

(iii) There are other sources of the nebular radiation - e.g. 
emitters from the hot star wind and ionizations by collisions. 
The former accounts for about 10% of the total EM for well 
developed winds (e.g. AG Peg) and the latter can be important 
for nebulae with T e <: 20 000 K. Therefore the emission mea- 
sure produced solely by emitters of the giant's wind is lower 
than what we measure. From this point of view our estimates 
of Mw represent upper limits as they were derived from the to- 
tal EM. For example, the mass-loss rate derived independently 
from accretion process in CQDra (Mw = 8.2 lO~ 8 M0yr~', 
Skopal 2005) is by factor of ~ 1 .5 lower than that derived here 
(M w = 1.2 10~ 7 M Q yr 1 , Table 3). 

Variability: The reconstructed SED at very different binary 
positions confirmed the orbitally-related variation in the emis- 
sion measure which causes the well-known wave-like variation 
in the LCs as a function of the orbital phase. The nature of this 
type of variability was recently investigated by Skopal (2001a). 
Here we used his conversion between the observed EM and the 
magnitude to reconstruct the U LCs of RW Hya and S Y Mus 
from the IUE spectra (Figs. 17 and 18). The observed ampli- 
tude of the wave-like variation is proportional to the ratio of 
fluxes from the nebula and the giant, which is a function of 
the wavelength - flux from the giant/nebula increases/decreses 
with increasing lambda. The contribution from the giant at 
the considered passband depends mainly on its spectral type, 
while the quantity of the nebular flux results from ioniza- 



tion/recombination processes and can be very different for in- 
dividual objects. Figure 25 shows that the contribution from 
the giant in the UBV bands is a function of its effective tem- 
perature, which determines the wavelength-dependence of the 
wave-like variation in LCs: Yellow objects show AU/AV » 1, 
while those containing late-type M-giants have AU/AV ~ 1 -3. 
In the column of the same r e tf , objects with a strong nebula are 
located at the bottom (i.e. AU/AV ~ 1). This is the case of 
V 1 329 Cyg at the column of r eff = 3 300 K and BF Cyg among 
objects at T eS = 3 400K. 

5.2.2. Stellar continuum radiation 

The hot star temperature, T^: It was not possible to deter- 
mine directly this parameter with our fitting procedure for 
Th ^ 10 5 K. As a result we selected appropriate values of Th 
according to the ionizing capacity of the hot source and the 
observed emission measure: 

(i) If the parameter 6 < 1 for T h ^ 10 5 K we adopted 
the Zanstra temperature for the ionizing object according to 
MNSV. This was the case of CI Cyg, SY Mus, and AG Peg. For 
two limiting cases, EG And( 10/09/95) and BFCyg( 10/04/92), 
we adopted the corresponding minimum temperatures, T™ n , 
because the referred Zanstra temperatures were below this 
lower limit. Note that for Th < T™ n the ionizing capacity is 
not sufficient to produce EM \, S - For RWHya we adopted Th 
= 110 000 K, because its ionized zone is very open even for 
Th = = 80 000 K and the suggested Zanstra temperature is 
below r™ n . In this special case it was also possible to get this 
temperature directly by fitting the well exposed spectrum from 
19/07/79. 

(ii) If the parameter 6 ~ 1 for a certain temperature Th > 
10 5 K, i.e. EM hs is very large, we adopted Th = T™ 11 . Then 
the parameters derived from the SED, Lh and Oh represent 
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their lower and upper limit, respectively (ZAnd, VI 329 Cyg, 
AG Dra, CQ Dra; Table 3). 

(iii) For two cases, V443 Her and AX Per, it was possible to 
determine the 7\ temperature directly as the fitting parameter of 
our procedure, because of its low value. However, a relatively 
strong nebular radiation in these systems yields the parameter 
8 » 1 . Therefore, in Table 3 we also give parameters for the 
case 8 = 1 (T h = T™ n ). 

The effective radius, R 1 ^: According to Eq. (6), is 
the radius of a sphere which has the same luminosity as the 
observed hot object seen under the angle 26^. According to 
present knowledge, accretors in symbiotic binaries are white 
dwarfs of a mean mass M WD ~ 0.6 Mq (Corradi et al. 2003). 
However, only in the case of CQ Dra, from our sample of the 
objects, is the /?jj ff radius comparable to that of a white dwarf 
(Table 3). Here the very small effective radius limits the white 
dwarf radius to Kwd ~ 0.006 Rq, which implies a rather mas- 
sive white dwarf, M W d ~ 1 Mq (Hamada & Salpeter 1961). 
Effective radii derived from the SED of other objects are a fac- 
tor of about 10 (or more) larger than a typical white dwarf's 
radius. Qualitatively, this can be understood as a result of accre- 
tion process that creates an optically thick photosphere around 
the mass core. Below we investigate this possibility: 

(i) The extreme cases of SEDs with 8 > 1 (i.e. T h < T™ n ) 
suggest a disk-like structure of the hot source. The observed 
hot stellar radiation is not capable of giving rise to the strong 
nebula in the system which constrains that an unseen part of the 
hot source radiates at a higher temperature to produce the nec- 
essary surplus of the ionizing photons. Thus, the inner ionizing 
source of the hot object has to be obscured by the accreted mat- 
ter more in directions forwards the observer than at its poles. 
The effect is significant for highly inclined systems. 

(ii) Models corresponding to 8 < 1 . If the hot stellar source 
radiation mimics that of a stellar photosphere, then it has to 



satisfy the LRT relation, which can be expressed in solar units 

as 

-(£hmSH'-(£)- 

A good agreement between this relation and our quantities of 
Lh and R 1 ^ from Table 3 confirms the abovementioned sugges- 
tion (Fig. 26). The best fit gives 41og(r h /r Q ) = 5.0 + 0.04, i.e. 
they radiate at a characteristic temperature of 103 000 + 2 500 K 
for Tq = 5 800 K. Here, to get a better solution, we omitted the 
cases of V1329Cyg and CQDra(l 1/10/86) (see Fig. 26). Thus 
we cannot observe just the accretor's surface. 

Variability: Observations at/around opposite conjunctions 
of the binary components confirmed the presence of an ad- 
ditional source of extinction, which attenuates the continuum 
(seemingly) at all wavelengths and is more pronounced at posi- 
tions around the inferior conjunction of the giant (here Z And, 
V1329Cyg). MNSV first noticed this extinction process. The 
effect is known for RWHya, SYMus, BFCyg (e.g. Dumm et 
al. 1999). Home et al. (1994) modeled the influence of the 
absorbing medium to the ultraviolet, optical and near-infrared 
spectrum of of the white dwarf in the eclipsing dwarf nova 
OY Car. They found that besides deep features around 1 600 
and 2 000 A the absorbing gas also produces a modest opti- 
cal depth in the Balmer and Pashen continuum which makes 
it lower by about 20% in the ultraviolet (see their Fig. 8). The 
magnitude of this effect is probably a function of the orbital 
inclination. We did not find it in the spectra of AG Dra and 
AG Peg. 
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5.3. The active phase 

5.3.1. The high-temperature nebula 

The HTN is seen directly during optical eclipses of all inves- 
tigated stars (BFCyg, CHCyg, CICyg, ARPav and AX Per) 
as the only component of the hot object radiation in the UV 
spectrum. Out of eclipses it is indicated by a non-zero level of 
the Rayleigh attenuated far-UV continuum. The nebular nature 
of this component was proved directly by the optical spectrum 
containing the Balmer jump in emission taken simultaneously 
with the ultraviolet measurements of the strong 1993-outburst 
of YYHer (Fig. 16). Also simultaneous UV/optical/radio ob- 
servations of BF Cyg taken at the totality of the optical eclipse 
during its 1989-92 active phase suggest the nebular origin of 
this component of radiation. In the UV it is given by its flat 
profile (Fig. 8), in the optical by the very negative colour in- 
dex, U - B ~ -0.7 (Skopal 1992), and in the radio the ob- 
served flux is increased to 0.84 mJy with respect to precedent 
values (Seaquist et al. 1993). The only counterpart to the radio 
emission is the HTN in the ultraviolet, which thus confirms its 
nebular nature. Its profile corresponds to high electron temper- 
atures, r e > 3 10 4 K (Table 4) and it is complementary to other 
components of radiation throughout the UV/optical/IR wave- 
lengths to match the observed SED. These characteristics imply 
that the HTN-radiation is not subject to the Rayleigh scattering 
process and the corresponding ionized medium is not eclipsed 
by the giant. Thus a major fraction of the HTN region is located 
at least ~ 1 R s above/below the orbital plane. The HTN appears 
to be very strong in the spectra from active phases. During tran- 
sitions to quiescence its emission measure gradually decreases 
(see figures for Z And, CHCyg, BFCyg). Other IUE observa- 
tions also support such evolution. For example, the HTN emis- 
sion measure from AX Per spectra on 30/11/90 (Fig. 21) de- 
creased by a factor of about 2.5 on its transition to quiescence 
on 04/09/94 (SWP52027 + LWP29094, <p = 0.01 - not pre- 
sented here). During quiescent phases it can be recognized as 
a very faint emission in UV (e.g. RW Hya, Fig. 17). Pereira et 
al. (1995) noted the presence of such the remaining light also 
during the eclipse of SYMus. They assigned it most likely to 
nebular emission. This suggests that creation of a strong HTN 
region is connected with the mass-outflow from the central star 
during outbursts. Below we discuss its origin. 

(i) Photoionization of a rather massive hot star wind can 
be excluded, because this would produce a significant emission 
within about 100 Rq from its source. The emissivity of such the 
wind is proportional to the square of the particle concentration, 
which is diluted in the wind with the radial distance r as 1/r 2 . 
In addition, a small wind velocity at the vicinity of its source 
makes the particle density even higher, which significantly en- 
hances the emissivity of the inner parts of the ionized wind (cf. 
Fig. CI of Skopal et al. 2002a). If this were the case, we would 
have observed a significant decrease in the emission from the 
HTN during the total eclipse. This, however, is not the case. 

(ii) Photoionization of the neutral particles from the giant's 
wind would produce nebular emission at T e ~ 1 - 2 x 10 4 K 
as we observe during quiescent phases. However, this does not 
correspond to the HTN properties. 
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(iii) The ejected material can make the surrounding neb- 
ular medium more opaque, which could lead to a significant 
heating of the outer layers of the symbiotic nebula, because 
the ionizing radiation is modified by absorption in the nebula - 
the higher-energy photons penetrate into the gas at larger (i.e. 
more opaque) distances from the central star, and thus the mean 
energy of the produced electrons is higher there (Osterbrock 
1974). The result probably depends strongly on the structure of 
the surrounding ionized material. 

(iv) In the analogy to the colliding wind model (e.g. Girard 
&Willson 1987; Nussbaumer& Walder 1993), a high temper- 
ature collisionally-ionized shock layer, which can be ionized 
also by the radiation field of the hot star, can be created well 
above the accretion disk. The collisions are between the high- 
velocity hot star wind and the slowly moving particles of the 
cool giant wind and/or circumbinary material. The temperature 
of this region can reach well above 10 6 K. Problems here are 
connected mainly with the efficiency of cooling the hot plasma 
and properties of the driven hot star wind. Nevertheless, phys- 
ical separation of the HTN from the central hot object, its very 
high electron temperature and production of the emission lines 
with very high ionization potential could be produced by this 
type of interaction. 

5.3.2. The low-temperature nebula 

The LTN is subject to eclipses, which determines its location in 
the vicinity of the central star and limits its linear size to 2 R g . 
This case is reproduced by CICyg, ARPav and AX Per, for 
which there are available ultraviolet observations taken in and 
out of their optical eclipses during activity. The radiative prop- 
erties of the LTN are characterized by its emission measure and 
luminosity. They are determined by the fitting parameters &n, 
T e and a in relations (13) and (14). Its electron temperature is 
between 10000 and 18 000K, which means that this region is 
predominately due to photoionization. Geometrical properties 
of the LTN - mainly its size and location - suggest its very high 
ionization degree. According to the ionization formula (e.g. 
Gurzadyan 1997), N+/N\ ~ a few x 10 3 even for very high 
electron concentrations of iV e ~ a few x 10 10 cirT 3 , because of 
a high value of the dilution factor (see below). Therefore and 
due to a small or negligible contribution from the doubly ion- 
ized helium (a ~ 0) we can assume N e = Nh+. Further, for 
the temperature of the ionizing source = 10 5 K, the av- 
erage electron temperature T e = 14000K, a dilution factor 
W = (R*/Rltn) 2 ~ (0.01/100) 2 = 10 8 and the optical depth of 
the nebula t c = k c Ni Run ~ 1 (* c = 6.3 10~ 18 cm 2 is the hy- 
drogen absorption coefficient at the Lyman limit, N\ is the con- 
centration of neutral H atoms and for /? LTN we adopt 100 Rq 
as a maximum), the ionization formula provides the electron 
density 

N e ~ 10 9 cnT 3 . (32) 

An independent estimate of 7V e can be made from the size of the 
LTN (it is subject to eclipse) and its emission measure. If we 
approximate its volume by a sphere of radius Rltn = 100/?q 
then its emitting volume V LTN = 1.4 10 39 cm 3 . Consequently, 
the emission measure obtained from the spectra out of eclipses, 
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EM LTN = N 2 e 

Vltn = 10 58 - 10 60 cm 3 , yields the average elec- 
tron concentration 

N e ~ a few x 10 9 cm" 3 , (33) 

which is comparable with the estimate made above using the 
ionization formula. This agreement supports that the LTN is 
mainly due to photoionization. Kenyon et al. (1991), based on 
an analysis of the line emitting regions in CICyg, also sug- 
gested the presence of two nebular regions in this system - a 
compact nebula located within the Roche lobe of the accretor 
and a low density, collisionally-ionized nebula having a form 
of bipolar streams above/below the disk's surface. 

5.3.3. The hot stellar source 

The region of the HSS is subject to eclipses (Figs. 8, 9, 10, 19 
and 21). Its radiation can be approximated by that of a black 
body at low temperatures. Below we summarize its radiative 
and geometrical properties as follows: 

1. Figure 26 shows: (i) The HSS radiation of most of ob- 
jects satisfies that of a stellar photosphere at a characteristic 
temperature of ~ 22 800 K. The HSS of CH Cyg and TX CVn 
radiates even at lower temperatures; we did not include them in 
the fit in Fig. 26. (ii) There is a significant enlargement of the 

radius by a factor of ~ 10 with respect to quantities from 
quiescence. This and the simultaneous presence of the two- 
temperature-type of the UV spectrum (see bellow, Sect. 5.3.4) 
suggest that the optically thick medium expands at the orbital 
plane in the form of a disk, (iii) The observed luminosity of the 
HSS during the lst-type of outbursts (see below; Fig. 26) de- 
creases by a factor of about 2-10 with respect to the hot object 
luminosity during quiescence (Table 3 and 4). This is because 
the HSS, which develops during these outbursts, blocks a frac- 
tion of radiation from the centre in the direction to the observer. 
We discuss this situation in Sect. 5.3.6 in more detail. 

2. The HSS radiation is attenuated by the Rayleigh scat- 
tering process for all the well-known eclipsing systems pre- 
sented here (BFCyg, CHCyg, CICyg, ARPav and AX Per; 
Table 4). The effect is measured at any orbital phase. We there- 
fore believe that the absorbing gas is located at the outer rim 
of the disk. On the other hand the presence of the Rayleigh- 
attenuated far-UV continuum in Z And and TX CVn confirms 
the recent suggestion of their high orbital inclination (Skopal 
2003a; Skopal et al. 2004), and thus suggests it also for CD- 
43° 14304 (Fig. 5). Generally, all the orbits of the objects show- 
ing the two-temperature-type of UV spectrum during activity 
should be highly inclined to the observer. Here this applies to 
AEAra(Fig. 4). 

3. Outbursts of AGDra and the major 1993-outburst of 
YYHer differ significantly from those described above. A 
strong nebular and very hot stellar components of radiation de- 
veloped immediately at the beginning of these outbursts. The 
hot stellar radiation is characterized by a higher temperature 
7h ~ 165 000 K, a larger Lh and smaller /?jj ff with respect to 
values from quiescence (Fig. 26). The case of AGDra could 
be caused mainly by an inclination effect: This system is seen 
more from its pole (i = 60 ± 8°, Schmid & Schild 1997) and 
thus the central ionizing source can be directly observed. As a 



result both a very hot and luminous stellar and strong nebular 
component of radiation can be measured directly. The observed 
large difference in the HSS temperature during the minor 1981 
and the major 1993 outburst of YYHer (Fig. 16) is probably 
related to different properties of the mass outflow. Material dur- 
ing the major outburst was ejected probably at very high veloc- 
ity, which precluded formation of a cool stellar pseudophoto- 
sphere. As a result the simultaneous massive injection of new 
particles (emitters) by the active star and the increase of its ion- 
izing capacity enhanced the nebular emission by a factor of 
-20 (Table 4). 

Accordingly, these observational properties of the HSS dur- 
ing active phases suggest two types of outbursts: The 1st type 
is characterized by a larger /^ ff and lower Th (points 1 and 2 
above) and the 2nd type is determined by a smaller R^ and 
higher 7h (the point 3 above) with respect to their quiescent val- 
ues (Fig. 26). All the active objects with a high orbital inclina- 
tion show characteristics of the 1st type of outbursts. Gonzalez- 
Riestra et al. (1999) identified 'cool' and 'hot' outbursts of 
AGDra. They lie within the 2nd type of outbursts. According 
to the point 1. (iii) above one can speculate that the lst-type of 
outbursts represent a special case of the 2nd-type due to a high 
orbital inclination (Sect. 5.3.6). Finally, we note that a large 
number of objects with highly inclined orbits, among those in- 
tensively investigated, reflects a systematic selection, because 
the symbiotic phenomenon is best recognizable for such sys- 
tems. 

5.3.4. The two-temperature UV spectrum 

The presence of the HTN, LTN and a warm HSS components 
of radiation determines the two-temperature-type of the ultra- 
violet spectrum. In our approach, the cool spectrum is pro- 
duced by the warm HSS radiation (~ 22 000 K) and the hot 
spectrum is represented by the emission from the HTN+LTN 
powered by an unseen hot ionizing source (<: 10 5 K) in the sys- 
tem. Superposition of these components results in the observed, 
more or less flat ultraviolet continuum. The far-UV region is 
dominated by the HSS radiation, while the near-UV continuum 
often contains a large contribution from the LTN. The slope of 
the HTN radiation throughout the ultraviolet is very small. This 
type of the U V spectrum represents a common feature of active 
phases classified here as the 1st type. 

5.3.5. Basic structure of the hot active object 

Radiative and geometrical properties of individual components 
of radiation observed during the 1st type of outbursts allow us 
to reconstruct a basic configuration of the hot active object. The 
following points are relevant. 

(i) The two-temperature-type of the UV spectrum suggests 
the presence of an optically thick disk-like torus around the hot 
star, which is seen approximately edge-on. The disk's outer rim 
has to be flared sufficiently to occult permanently the central 
source of the hot ionizing radiation, which is not seen directly 
in the spectrum. Under such conditions, the energetic radia- 
tion from the central source is absorbed and diffused by the 
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Fig. 27. Schematic representation of a basic structure of hot objects during outbursts of the 1 st-type. The left panel shows example of individual 
components of radiation we isolated from the ARPav spectrum observed on 10/05/81. Note that the HSS radiation is Rayleigh attenuated at the 
far-UV region and veiled by a forest of blended Fe n absorptions. The right panel shows a sketch of the corresponding emitting regions as seen 
on a cut perpendicular to the orbital plane containing the accretor. During outbursts the disk's outer rim has to be flared to occult permanently 
the central source of the hot ionizing radiation. The model is described in Sect. 5.3.5. 



disk body in the direction to the observer. At the distance Russ 
from the center, it releases the disk at much lower tempera- 
ture (~22000K). We identify this warmer pseudophotosphere 
as the HSS. The HSS radiation then continues throughout the 
neutral gas surrounding the accretor beyond the disk, which 
attenuates its amount by the Rayleigh scattering process and 
affects (often drastically) its profile by numerous absorptions 
of mainly Fen lines, the so-called 'iron curtain'. 

(ii) The geometrical and physical properties of the LTN lo- 
calize this region above/below the disk plane approximately 
within the linear size of the HSS (it is also subject to eclipse). 
It is probably also flared at the outer edge of the disk, because 
directions from the hot inner parts of the disk to poles have to 
be rather free for high-velocity particles and ionizing photons, 
which give rise the uneclipsed emission from the HTN. In ad- 
dition, the LTN has to 'see' well the ionizing radiation from 
the hot center. This view is supported by theoretical modeling 
of Plavec & Hubeny (1994), who found that the vertical outer 
edge of accretion disks in symbiotic binaries can be extremely 
extended. In the real case the vertical extension is limited by the 
eclipsing nature of the LTN. For the sake of simplicity, we draw 
schematically the LTN vertical dimension as a linear function 
of the disk radius (cf. Fig. 27). 

(iii) The HTN region has to be located high enough 
above/below the disk (> 1 R g ) not to be eclipsed by the gi- 
ant (e.g. Fig. 10 and 19). However, our observations were not 
sufficient to determine its geometry more accurately. The HTN 
can be in part collisionaly ionized, which suggests that a mech- 
anisms similar to that of the colliding wind model could be 
responsible for its rise (Sect. 5.3.1). 

A sketch of the above described geometrical structure of the 
hot object in symbiotic binaries during the 1 st-type of outbursts 
is shown in Fig. 27. 



5.3.6. Luminosity of the hot active object 

In this section we discuss the total luminosity of hot objects 
undergoing the 1 st-type of outbursts. The edge-on disk/torus, 
which develops during these outbursts, blocks a fraction of ra- 
diation from the accretor in the direction to the observer. We 
measure directly only its part in the form of the HSS radia- 
tion. The remaining part of the accretor's luminosity can be- 
come visible due to its interaction with the medium surround- 
ing the hot object. We detect only its fraction converted into the 
nebular emission in the form of the LTN and HTN radiation. 
However, a part of the accretor's luminosity, L out , can escape 
the system. Therefore we write the bolometric luminosity, L T , 
produced at/around the white dwarf as 

Lj = LuSS + ^HTN + ^LTN + ^out- (34) 

The not-detectable light includes the radiation, which is not ca- 
pable of ionizing the surrounding medium (for Th = 10 5 K and 
hydrogen it is about 10% of the total radiation (Fig. 2 of Skopal 
2003b)), and a part of ionizing photons, which are not con- 
verted into the nebular emission. Thus we, generally, observe 

< s « Li. (35) 

This is the case of all systems showing signatures of the 1st type 
of outbursts (compare data from Tables 3, 4 and Fig. 26). On 
the other hand when viewing the system more pole-on we can 
detect directly all the radiation released during outbursts. An 
example here is AGDra. During its 'hot' outbursts (Sect. 4.13) 
we have 

L° bs * L T = Iq, (36) 

if one adopts the total luminosity as that from quiescence, Lq. 
Note that for Th > T™" 1 one can easily get equivalence be- 
tween the luminosities in the above relation. This suggests 
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an accretion-powered nature of these outbursts. However, for 
'cool' outbursts we measure 

K hS » £q. (37) 

which signals an additional source of the energy during these 
events, probably from a TNR. Some details on the nova-like 
eruptions of symbiotic stars can be found in Mikolajewska & 
Kenyon (1992b). 

5.4. Comparison with previous models 

A modeling of the composite continuum for a large sample of 
well-studied symbiotic stars similar to the present one has been 
carried out by MNSV and KW. The basic conceptional differ- 
ences between these and our approach were introduced in the 
section 1. Here we present a concise summary of differences 
between our models and those suggested by other authors. The 
following points are relevant: 

(i) Our SED model (Eq. 13) is less sophisticated than that 
suggested by MNSV. We did not include a geometry of radi- 
ating regions and consider just basic sources of radiation and 
effects that can be recognized by observations. This simplifies 
significantly modeling the SED and makes it more accessible 
for a wider community. On the other hand, this does not result 
in decisive differencies in and Lh from quiescent phases. 

(ii) In our approach the electron temperature represents the 
fitting parameter, while MNSV considered it to be within a 
small range (Sect. 1). This revealed drastic variations in T e even 
for the same object (e.g. ZAnd, AG Peg). This allowed us to 
determine corresponding EMs more accurately. We note that 
an adoption of a 'typical' value of T e can result in an incorrect 
SED (e.g. Fig. 4 in Munari et al. 1997a). 

(iii) In our approach we used temperatures required by the 
basic photoionization model (i.e. T™ 11 (Eq. 23) or r Zanstra , see 
Sect. 5.2.2) as nebulae during quiescent phases are largely ra- 
diatively ionized (r e :S 20 000 K). In cases of AX Per and 
V443 Her a large temperature difference, Ar = T™ n - T m ~ 
35 000 K (Table 3), suggests the presence of a disk-like struc- 
tured material around their accretors even during quiescent 
phases (Sects. 4.15 and 4.21). MNSV discarded the presence 
of accretion disks in their solutions, whereas KW considered 
them as a reasonable possibility for some cases. Most recently 
Sion et al. (2002), Sion (2003) and Kolb et al. (2004) ap- 
plied model atmospheres in spectral analyses of the hot stellar 
component of radiation from RW Hya, AG Dra and EG And. 
They also indicated significantly lower effective temperatures 
than those given by conditions of photoionization in these sys- 
tems. Effective temperatures of their models and our T™ n for 
EG And and RW Hya (Table 3) give a similar Ar to what we 
found for V443 Her and AX Per. This suggests that we detect 
temperatures of different regions in the hot object. Within our 
interpretation we measure directly a cooler pseudophotosphere 
(~ edge-on disk-like structured accreting material), while the 
underlying hotter accretor's surface can be 'seen' only by the 
surrounding medium, which converts its radiation to the neb- 
ular emission, based on whose properties we can derive indi- 
rectly the temperature of the hottest part of the ionizing source 
as r Zanstra and/or T™ n . Such a structure of hot components can 



explain the apparent inconsistence between effective tempera- 
tures from model atmospheres analyses and those required by 
the observed nebular emission. In accordance with this view 
the temperature difference Ar should be larger for objects with 
a high orbital inclination and large mass-loss rates, and vice 
versa. This idea could be tested by applying the model atmo- 
sphere spectral analysis to objects seen rather pole-on and with 
a small Mw (e-g- CQ Dra in our sample). Such testing is beyond 
the possibilities of our simple black-body fitting. 

(iv) The mass-loss rate was calculated for a more realis- 
tic structure of the giant's wind. We considered a /Haw wind, 
whereas MNSV adopted a constant wind. In addition, we used 
a more accurately dete mined EMs (the point (ii) above). This 
difference, however, did not yield a large quantitative effect. 
Mikolajewska et al. (2002b) estimated M w from their radio 
mm/sub-mm observations assuming a partially ionized red gi- 
ant wind determined by opened ionized zones. This approach 
is, in principle, similar to our calculations. A comparison of 
our M w quantities with those mentioned above is satisfactory 
(Fig. 24). 

(v) We quantified contributions from cool components by 
comparing an appropriate synthetic spectrum to IR flux-points 
given by photometric measurements. KW reconstructed the 
spectrum of cool components according to energy distributions 
of Johnson (1966) and Lee (1970) and scaled it to absolute V 
magnitudes of K3, M2 and M4 giants. MNSV restricted their 
analysis only to UV spectra. Our approach allowed us to esti- 
mate effective temperatures and angular radii of giants indepen- 
dently from their spectral types. Application of synthetic spec- 
tra of red/yellow giants is of a particular importance to model 
SEDs for objects with a considerable cool component contribu- 
tion at the near-UV (yellow symbiotics and those with a small 
nebular emission, e.g. CQ Dra) to get a correct T e and EM. For 
comparison, a simple blackbody fit gives a large contribution 
at these wavelengths (e.g. Fig. 7 in Skopal et al. 2000). 

(vi) A major qualitative difference is that we modeled spec- 
tra from all the quiescent phases, active phases and eclipses, 
whereas MNSV restricted their analysis exclusively to quies- 
cent phases at binary positions well outside eclipses. Their tri- 
als to model active phases failed (Fig. 8 and Table 6 in MNSV). 
Application of our method to spectra taken during active phases 
revealed their common properties (Sect. 5.3; Fig. 27). This rep- 
resents the main novelty of the present paper. For compari- 
son, in our model of the active hot object of the lst-type the 
flared disk permanently occults its hotter inner parts in the 
line-of-sight, whereas according to previous suggestions (e.g. 
KW; Kenyon et al. 1991; Mikolajewska & Kenyon 1992a) 
we observe radiation from the inclined disk's surface including 
its central hot parts. However, fits to the observed continuum 
by using parameters of the latter models are poor (e.g. Fig. 9 
in MNSV). A quantitative application of standard disk mod- 
els to the UV/optical/IR continuum was carried out by Skopal 
(2003b) for ARPav. This approach failed to explain the UV 
continuum with the disk and its boundary layer. 

(vii) To select the continuum flux-points we considered 
the influence of the iron curtain. This led to determination of 
more realistic model parameters for spectra from mainly active 
phases. Particularly, we obtained solution also for CHCyg, on 
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which approaches of both the groups (KW and MNSV) failed. 
A first suggestion was made by Shore & Aufdenberg (1993) 
for the recurrent nova TCrB. However, ignoring the effect of 
the iron curtain for this case, Stanishev et al. (2004) obtained 
an ineffectual SED corresponding to false parameters, which 
are not consistent with the observed hot object luminosity. Also 
Skopal (2003a) obtained false quantities of He ++ abundance 
and T e of the LTN in his modeling the Z And spectra from the 
active phase. 

(viii) Finally, the presented comprehensive atlas of mod- 
eled and observed SEDs for intensively studied S-type sym- 
biotic stars accompanied by LCs represents a large source of 
information for future works. 

6. Conclusions 

We introduced a method of disentangling the composite spec- 
trum consisting of two stellar components of radiation (from 
cool giants and hot objects in the binary) and nebular radia- 
tion. We applied the method to 21 S-type symbiotic binaries 
to isolate individual radiative components from their 0.12 - 
5 /J.m SED in the continuum during quiescent as well as ac- 
tive phases. The modeled continuum is determined by 8 fitting 
parameters: Temperatures, 7^, 7h and T e , angular radii, 6 g and 
#h, scaling of the nebular emission, fc N , hydrogen column den- 
sity, «h and average abundance of doubly ionized helium, a 
(Sect. 3.1). These parameters then determine luminosities and 
radii of the stellar components and the emission measure of the 
nebula. Common properties of individual radiative components 
and their sources can be summarized as follows: 

1. Radiation from the giant: (i) The SED of giants is 
defined by the fitting parameters, T s g and 6 g . Those having 
independently-determined radii (11 objects from our sample, 
Fig. 23) then provide the distance from their 6 g . Relations (27) 
and (28) show that their characteristics satisfy best those of nor- 
mal giants. This confirms a generally accepted view that cool 
components in S-type systems are normal giants as given by 
their spectral types (e.g. Miirset & Schmid 1999). (ii) The ef- 
fective temperature of giants is critical for the amplitudes of the 
wave-like orbitally-related variation in the LCs (Fig. 25). 

2. Hot object radiation during quiescence: (i) The nebu- 
lae radiate at a mean temperature of 19 000 K. By compar- 
ing the measured and calculated emission measure we derived 
the mass-loss rate from giants via their wind as Mw = a few 
x 10" 7 M Q yr 1 (Table 3, Fig. 24). (ii) Stellar continuum radi- 
ation determines the effective radii of the hot objects to be a 
factor of ~ 10 larger than the accretors. We thus observe pseu- 
dophotospheres due to the accretion process (Fig. 26). The ex- 
treme cases, for which the SED suggests T h < T™ n (the ob- 
served hot source is not capable of giving rise to the nebular 
emission), signal the presence of a disk-like structure of the 
circumstellar matter around the accretor. This is in agreement 
with theoretical modeling of accretion from the giant's wind 
that predicts formation of a stable, permanent disk around the 
accretor (Mastrodemos & Morris 1998). 

3. Hot object radiation during activity: (i) The stellar radi- 
ation satisfies that of a black-body photosphere. There are two 
types of outbursts: The lst-type with the HSS radiating at sig- 



nificantly lower temperatures than during quiescence and the 
2nd-type with higher characteristic temperatures (cf. Fig. 26). 
All the active objects with a high orbital inclination show fea- 
tures of the 1st type of outbursts: Rayleigh attenuation at any 
orbital phase, a strong influence from the iron curtain absorp- 
tions, a low characteristic temperature of ~ 22 000 K, strong 
HTN emission, two-temperature-type UV spectrum and en- 
largement of R 1 ^ by a factor of ~ 10 with respect to the qui- 
escent values. The last two characteristics suggest an expan- 
sion of an optically thick medium at the orbital plane in the 
form of a disk, (ii) The low-temperature nebula (LTN) radi- 
ates at a mean electron temperature of 14 000 K. This region is 
subject to eclipses, which restricts its maximum linear radius 
to 2R g . Its emission measure, 10 58 - 10 60 crrT 3 , suggests the 
electron concentration /V e ~ a few x 10 9 ctrT 3 . (iii) The high- 
temperature nebula (HTN) is characterized by T e > 30 000 K 
and its major fraction is located well away from the hot ob- 
ject, not subject to eclipses. Its emission measure is between 
10 58 - 10 61 cm . Its origin could be within a scenario of the 
colliding wind model. 

The characteristics of the UV-continuum and its compo- 
nents we observe during the 1st type of outbursts can be in- 
terpreted in terms of an edge-on flared disk surrounded by the 
neutral material at the orbital plane and the nebulae located 
above/below its surface. The model is described in Sect. 5.3.5 
and in Fig. 27. 
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